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SUMMARY 


The  effects  of  rain  and  fog  are  examined  for  three  transmissometer 
and  three  radar  systems:  double-ended  (direct  path)  with  square-wave 
modulation  and  square-law  detector,  single-ended  (folded  path)  with 
square-wave  modulation  and  square-law  detection,  single-ended  CW  with 
quadrature  receiver,  single-ended  pulsed  with  post -detect ion  pulse 
integration,  single-ended  pulsed  with  quadrature  receiver,  and  single- 
ended  with  linear  FM  pulse  compression.  Computer  algorithms  were 
developed  to  find  the  probability  of  detection  as  a  function  of  range 
for  the  radar  systems,  and  the  signal -to-noise  ratio  as  a  function  of 
range  for  the  transmissometer  systems.  These  algorithms  are  used 
recursively  to  determine  the  range  at  which  specified  performance  is 
attained  as  a  function  of  rain  rate  or  fog  liquid  water  content. 
Calculated  results  are  presented  for  exponential  drop-size 
distributions;  provision  has  been  made  for  the  extension  to  other 
distributions. 
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SECTION  I 


INTRODUCTION 

The  basic  purpose  of  this  study  is  to  investigate  the  performance 
of  transmissometers  which  could  be  used  to  measure  the  transmission  over 
a  path  on  which  millimeter  wave  radars  are  simultaneously  evaluated. 

For  the  most  part,  a  detailed  analysis  of  six  different  types  of  trans¬ 
mission  systems  was  performed  for  frequencies  of  35,  94,  140,  and  240 
Gigahertz  (GHz).  Three  of  these  systems  were  transmissometers  and  three 
were  radars  that  might  be  evaluated  on  the  range.  These  frequencies 
have  been  selected  because  they  represent  the  "windows"  of  low  molecular 
absorption  as  shown  below  by  Figure  1. 

In  this  project  two  major  problems  were  examined.  First,  we  con¬ 
sider  the  problem  of  accurately  modeling  the  attenuation  through  rain 
and  fog  due  to  scattering  and  absorption.  These  effects  are  character¬ 
ized  in  terms  of  frequency  and  meteorological  quantities  such  as  rain 
and  fog  drop-size  distributions,  rain  rate,  temperature,  and  liquid 
water  content. 

Second,  the  six  systems  were  analyzed  in  terms  of  the  common  per¬ 
formance  parameters,  false  alarm  probability  and  probability  of  detec¬ 
tion  for  the  radars  and  signal -to-noise  ratio  for  transmissometers. 
Equations  which  predict  the  performance  parameters  at  a  specified  range 
are  derived.  Then  these  equations  are  used  to  find  the  effect  of 
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Figure  1.  Specific  Attenuation  In  clear  air  versus  frequency, 
(after  Rosenblum  [1].) 
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weather  on  this  range  at  the  various  frequencies. 

The  term  "double-ended"  will  be  used  to  describe  a  system  in  which 
transmitter  and  receiver  are  separated  and  direct  transmission  occurs 
between  them.  A  "single-ended"  system  is  one  in  which  the  signal  Is 
reflected  by  a  target  to  the  receiver,  which  is  located  close  to  the 
transmitter. 

In  this  project,  the  following  transmissometer  and  radar  systems 
were  examined: 

(1.)  Double-ended  transmission  path  with  audio  (square-wave) 
modulation  and  square-law  detection.  Note  this  cannot  be  used 
as  a  radar.  Analysis  is  not  sensible  for  probability  of 
detection  because  there  is  no  target.  Instead,  a  detailed 
analysis  in  terms  of  signal-to-nolse  ratio  is  performed. 

(2.)  Single-ended  transmission  path  with  audio  (square-wave) 
modulation,  and  square-law  detection  using  a  calibrated  tar¬ 
get.  This  is  a  transmissometer  which  could  be  used  as  a 
radar  if  range  information  is  not  required.  Analysis  in  terms 
of  signal -to-noise  or  probability  of  detection  could  be  per¬ 
formed,  however,  signal -to-noise  ratio  analysis  is  more 
meaningful . 

(3.)  Single-ended  transmission  path  with  audio  modulation, 
quadrature  detection  and  a  calibrated  target.  Similar  to  the 
previous  system,  this  is  a  transmissometer  which  could  be  used 
as  a  radar  if  range  information  is  not  required.  Once  again, 
signal -to-noise  ratio  analysis  techniques  are  more  useful, 
however,  a  probability  of  detection  and  probability  of  false 
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alarm  can  be  defined  since  we  do  have  a  calibrated  target. 

(4.)  Post -detect ion  (boxcar)  integration  over  a  single-ended 
transmission  path.  This  is  a  true  radar.  It  could  be  used  as 
a  relatively  inefficient  transmissometer  by  ranging  on  a  cal¬ 
ibrated  target.  For  this  system,  analysis  in  terms  of  proba¬ 
bility  of  detection  and  probability  of  false  alarm  is  useful. 
(5.)  Quadrature-detection  over  a  single-ended  transmission 
path.  As  with  the  previous  system,  this  is  a  true  pulsed 
radar  which  could  be  used  as  a  relatively  inefficient  trans¬ 
missometer.  Once  again,  probability  of  detection  and  proba¬ 
bility  of  false  alarm  is  most  useful  for  analysis  of  this 
system. 

(6.)  Pulse  compression  with  quadrature  detection  over  a 
single-ended  transmission  path.  This  radar  system  is  Identi¬ 
cal  to  the  previous  system  except  pulse  compression  Is  used. 
This  system  can  also  be  used  as  a  transmissometer.  As  with 
the  last  system,  probability  of  detection  and  probability  of 
false  alarm  is  most  useful  for  analyzing  this  system. 

The  basic  format  of  this  thesis  is  such  that  expressions  for 
signal -to-noise  ratio  are  given  in  Section  II.  A  derivation  of  expres¬ 
sions  relating  probability  of  detection  to  signal -to-noise  ratio  is 
given  in  Section  III.  Fquations  for  calculating  drop-size  distributions 
and  for  determining  the  attenuation  due  to  rain  and  fog  are  derived  in 
Section  IV.  The  results  obtained  from  computer  models  using  the  equa¬ 
tions  derived  in  Sections  II  to  TV  are  described  in  Section  V. 
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SECTION  II 


DERIVATION  OF  SIGNAL -TO -NOISE  EXPRESSIONS 

In  this  section,  the  signal -to-noise  ratio  expressions  for  the  six 
configurations  of  instrumentation  discussed  in  Section  I  are  derived. 

In  Section  II. A,  the  signal -to-noise  ratio  expressions  are  derived  for  a 
point  in  the  receiver  prior  to  any  demodulation  or  nonlinear  signal 
processing.  Therefore,  the  expressions  for  signal -to-noise  ratio 
correspond  to  the  output  signal  of  the  intermediate  frequency  amplifier. 
These  signal -to-noise  expressions  depend  upon  whether  a  one-way  or  a 
two-way  transmission  path  is  considered,  however,  they  are  independent 
of  the  type  of  signal  detection.  The  detection  process  is  considered  in 
Section  II. B.  The  expressions  given  in  Sections  II. A  and  II.B.l  to 
II.  B.5  are  based  on  Report  AFWAL-TR -81-1281  [2], 

A.  Deri  vat ion  of  Pre -Detect ion  Signal -to-Noise  Ratio  Expressions. 

Each  transmission  configuration  studied  may  be  described  in  terms 
of  one  of  the  following  groups: 

(1.)  The  system  uses  signal  transmission  over  a  one-way  path. 

(2.)  The  receiver  Is  located  next  to  the  transmitter  and  the 
signal  from  the  transmitter  is  reflected  from  a  calibra¬ 
ted  target  back  to  the  receiver.  This  is  a  radar  system 
but  possibly  without  a  requirement  for  range  resolution 
(transml ssometers . ) 

Let's  consider  the  derivation  of  the  signal -to-noise  expressions  for 


each  transmission  system  separately 


For  a  one-way  path,  the  power  density  at  the  receiving  antenna  is 


Ve= 


(1) 


where  P^  *  transmitted  power, 

Gt  =  transmitter  antenna  gain  in  direction  of  the  receiver, 
and 

R  =  distance  between  transmitting  and  receiving  antennas. 


Therefore,  the  power  received  is  given  by 

P, 


.Sits** 


(2) 


where  Pr  =  power  received  by  receiving  antenna, 
and  Ae  *  effective  aperture  of  receiving  antenna. 


The  effective  aperture  of  the  receiving  antenna  is  in  units  of  area  and 
may  be  defined  as 


Ae 


Gr 

*  “£r 


(3) 


where  Gr  =  gain  of  the  receiver  antenna  in  the  direction  of  the  trans 
mitter, 

X  =  wavelength  of  the  signal. 


The  Friis  transmission  equation  Is  obtained  by  combining  Equations  (2) 
and  (3). 


Pr 


(♦> 


If  losses  in  the  transmitter  and  receiver  system,  as  well  as 


attenuation  factors  such  as  rain  or  fog  in  the  transmission  path,  are 


Included,  then  the  received  power  Is  given  by 

-c* 


where  a  =  specific  attenuation  in  decibels  per  unit  length, 

Lt  =  transmitter-to-transmittlng  antenna  loss  (expressed  as 
efficiency). 
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Lr  =  receiving  antenna -to-recei ver  loss  (expressed  as  efficien¬ 
cy.) 

Now,  if  the  thermal  noise  power  at  the  pre-detection  filter  output 

signal  (usually  the  intermediate  frequency  (IF)  amplifier)  is 

*  =  kTpRF  (6) 

where  k  =  Boltzman's  constant  (1.38x10-23  Joules/°K), 

To  =  reference  noise  temperature  (290  °K), 

8  =  effective  pre-detection  receiver  filter  bandwidth  (Hz), 

F  =  receiver  noise  figure  referenced  to  290  °K, 

an  expression  for  the  signal -to-noise  ratio  in  the  output  signal  of  the 

intermediate  frequency  can  be  found.  This  signal -to-noise  ratio  for  a 

one-way  transmission  path  is  determined  by  combining  Equations  (5)  and 

(6).  Thus,  for  propagation  over  a  one-way  path,  the  signal -to-noise  is 

-aR 

ITT 

S  =  PtGtGrX2LtLr  10  .  (7) 

N  "T45)2JmcTofeF 

Now,  consider  systems  of  the  radar  type  transmitting  over  a  two-way 
path.  The  power  density  at  the  target  is  given  as 

-oR 

“ITT 

St  =  Mf  10  .  (8) 

TttRZ 

If  the  radar  cross  section  of  the  target  is  defined  as  a.  In  units  of 
area,  the  power  reradiated  in  the  receiver  direction  is 


-aR 

“TIT 

P'  =  PtGto  10  .  (9) 

TrrR7- 


Since  propagation  is  over  a  two-way  path,  the  power  density  of  the  echo 


signal  at  the  receiving  antenna  Is 

-2oR 

“ITT 


(10) 


The  received  power  is  obtained  by  multiplying  Sr  by  the  effective 


aperture  area.  By  use  of  Equation  (3),  the  radar  equation  becomes 

-2  aft 

“ITT 

Pp  =  PfG+GPaX2  10  .  (11) 

Including  effects  of  the  losses  in  the  transmitting  and  receiving  system 


-2cft 

IT 

Pr  =  PtGtGPLfLrX2g  10  .  (12) 

(4ir)JK4 

Therefore,  from  Equations  (6)  and  (12),  the  pre-detection  slgnal-to- 
noise  ratio  for  the  radar  transmissometer  is 

-2cft 

~I7T 

S  =  PtGtGPLtLroX2  10  .  (13) 

N  ■Tf«)3R*RTo'6T 

Equation  (7)  gives  the  pre-detection  signal -to-noise  ratio  for  a 
one-way  transmission  path,  whereas  Equation  (13)  is  for  a  two-way 
transmitter-receiver  configuration. 

B.  Derivation  of  the  Processed  Signal -to-Noise  Ratio  Expressions. 

The  six  configurations  described  In  Section  I  are  examined  to 
determine  the  signal -to-noise  ratio  after  detection  and  signal 
processing. 

1.  One-way  Transmission  and  Square-Law  Detection. 

This  configuration  transmits  and  receives  a  chopped  signal  over  a 

one-way  propagation  path.  As  shown  by  the  block  diagram  of  the  receiver 

In  Figure  2,  the  received  signal  Is  heterodyned  to  form  an  intermediate 

frequency  signal,  amplified,  detected  by  a  square-law  device,  and 
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filtered  by  a  narrow-band  (audio)  filter.  The  signal -to-noise  ratio 
expression  given  by  Equation  (7)  applies  to  r(t). 

The  chopped  signal  may  be  approximated  as  an  amplitude-modulatea 

signal  with  a  modulation  index  (u)  of  100%  since  only  the  fundamental 

harmonic  component  is  selected  by  the  audio  filter.  Therefore,  the 

input  signal  to  the  square-law  device  is  given  by 

r(t)  =  Ac[l+ycoswat]cos<DCt  +  n(t)  (14) 

where  Ac  =  carrier  amplitude, 
u  =  modulation  index, 

o)c  =  intermediate  amplifier  angular  frequency, 
ua  =  angular  center  frequency  of  audio  filter. 

The  factor  n(t)  represents  narrowband  Gaussian  noise  with  a 
spectral  density  determined  by  the  spectral  response  of  the  intermediate 
frequency  amplifier  and  centered  at  the  intermediate  angular  frequency 
“c*  Typical  values  for  wa/2ir  and  are  1  KHz  and  1  GHz 

respectively. 

With  an  input  signal -to-noise  as  given  by  Equation  (  7),  the 

signal -to-noise  power  ratio  at  the  audio  filter  output  signal  Is 

(E/N)  =  u2CS/Nl2[81f/Ba]{l+CS/N^(2+M2)}-1  (15) 

where  S/N  =  pre-detection  signal -to-noise  ratio  given  by  Equation  (7), 
Bif  =  equivalent  rectangular  bandwidth  of  the  intermediate  freq¬ 
uency  amplifier,  same  as  B  in  Equation  (6), 

Ba  =  bandwidth  of  the  audio  filter. 

Finally,  the  expression  for  signal -to-noise  ratio  of  a  system  with 
a  chopped  waveform,  modulation  index  of  100%,  and  square-law  detection 
is 

(E/N)  =  (16) 

where  (S/N)  is  given  by  Equation  (7). 
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2.  Two-way  Transmission  and  Square-Law  Detection. 

For  this  system,  the  receiver  is  identical  to  that  of  the  system 
with  one-way  transmission  and  square-law  detection  previously  discussed. 
However,  in  this  case  a  reflective  target  produces  the  received  signal. 
The  processed  signal-to-nolse  ratio  is  again  given  by  Equation  (16), 
but  with  (S/N)  defined  by  Equation  (13). 

3.  CM  Signal  Transmission  with  a  Quadrature  Receiver. 

This  receiving  configuration  consists  of  the  transmission  of  * 

CM  signal  over  a  two-way  path  with  the  use  of  a  reflecting  target. 
However,  a  quadrature  detector  is  used  rather  than  a  square-law  device 
as  in  the  previous  system.  A  block  diagram  for  a  quadrature  receiver 
is  shown  in  Figure  3,  where  u(t)  is  the  unit  step  function.  This  unit 
step  function  is  Included  only  to  show  that  all  integrators  are  reset  at 
t=0  (i.e.,all  earlier  signals  are  disregarded.) 

If  the  integrators,  which  are  not  continuously  operating  devices, 
are  replaced  with  equivalent  continuously  operating  filters  as  discussed 
below,  then  the  unit  step  function  u(t)  should  be  omitted.  In  this 
case,  the  detector  input  signals  are  identical  for  Figures  2  and  3.  For 
mathematical  convenience,  the  integrators  may  be  used  when  deriving  the 
output  signal -to-noise  ratio. 

Following  Buyukdura  [2]  and  Schwartz  [3]  to  obtain  the  output 
signal -to-noise  ratio,  a  single  correlation  detector  is  considered  as 
shown  in  Figure  4a.  The  Integrator  output  signal  is  given  by 

T 

y (T )  =  /r(t)s(t)dt  . 

0 


(17) 


Now,  consider  the  matched  filter  shown  In  Figure  4b.  The  filter  output 
signal  at  t=T  is  given  by  the  convolution 

y (T)  =  r(t)u(t)*h(t),  (18) 


that  is, 

y(T)  «  /r(t)u(t)h(T-t)dt.  (19) 

-00 


However,  the  definition  of  the  unit  step  function  and  the  causality 
requirement  for  the  filter  requires 

u (t )  =  0  ,  t<0  (20) 

and 


h(T-t)  =  0  ,  t>T  .  (21) 

The  solution  of  Equations  (20)  and  (21)  with  Equation  (19)  gives 

T 

y(T)  =  /r(t)h(T-t)dt  .  (22) 

0 


Note  that  if  the  filter  response  satisfies 

h(T-t)  =  s (t )  ,  0<t<T  (23) 


or  that 


h(t)  =  s(T-t)  ,  t>0  (24) 

then  y(T)  is  the  same  as  for  the  correlation  detector  of  Figure  4a. 
From  this  figure,  let's  consider  that  the  signals  are 

r(t)  =  Acoswct  (25) 

and 


s(t)  =  Ccosu)Ct  .  (26) 

From  Equation  (22),  the  filter  impulse  response  which  makes  it 
equivalent  to  the  integrator  network  must  then  be 

h(t)  =  Ccoswc(T-t)  .  (27) 
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If  T  is  chosen  to  be  an  integer  multiple  of  the  period  of  r(t), 
then  Equation  (26)  is  the  impulse  response  of  a  matched  filter.  That 
is,  when 

T  =  2nn  n  =  0,±1,±2,.  .  . 

03c 

then  Equation  (26)  becomes 

h(t)=  Ccos(2nn-u>ct)  =  Ccosoict  =  s(t).  (28) 

Imposing  the  restrictions  placed  by  Equations  (20)  and  (21),  one 
obtains 

h(t)  =  Ccos(»ict  u(t)  u(T-t)  .  (29) 

To  specify  the  transfer  function  of  the  filter,  we  must  find  the  Fourier 
transform  of  the  impulse  response.  This  can  be  found  from 

T 

H(ju)  =  C/cos  ucte-i^dt.  (30) 

0 

Applying  Euler's  Identity  for  cosu>ct, 

T 

H(joj)  =  C  /[ej^ct+e^J^ctle-i^dt 
7  0 

or 

T 

H(j<o)  =  C/[ej(“+“c)t+ej(u>-“c)tldt. 

7  0 

Performing  the  integration  shows  the  transfer  function  of  the  filter  to 
be  of  the  form 

j(oj+o)c)'7  T  j(<i>-(oc)7  T 

H(juj)  =  CT  e  sin(uH-a)c)?  +  e  sin(h>-(oc)?  .(32) 

(o.-a,c)T 

Because  the  matched  filter  shown  in  Figure  4b  is  linear,  the  signal 
and  noise  components  of  the  output  signal  may  be  examined  separately. 


The  signal  at  the  output  can  be  obtained  directly  from  the  correlation 
detector  of  Figure  4a.  and  Equations  (25)  and  (26)  as 

T 

Y(T)  =  /r(t)s(t)dt  (33) 

0 

T 

=  /ACcos2wctdt 

n 

Y(T)  =  ACT  .  (34) 

2 

This  expression  for  Y(T)  can  he  used  to  define  an  expression  for  the 
signal  power  (normalized  to  a  characteristic  Impedance  of  1  ohm)  by 

Ps  =  [ACT/2]2  .  (35) 

According  to  Haykin  [4]  and  Taub  and  Schilling  [5],  the  noise  power 
at  the  output  is 

00 

pny  =  1  /Sny(u>)du>  (36) 

2  7T  -oo 

where  Sny(a))  is  the  power  spectral  density  of  the  noise  component  of 
Y(T).  It  is  given  in  terms  of  the  power  spectral  density  of  the  noise 
component  at  the  input  by 

Sny(<o)  =  Snr|H(u>)|2.  (37) 

If  one  assumes  a  filter  with  a  rectangular  bandwidth  B  of  unit  ampli¬ 
tude,  the  noise  power  referred  to  the  output  is  given  by 

oo 

pny  =  _L_  /  $nr ! ^rect (U))f ^  da>  =  kTfjBF  (38) 

2  it  -» 

Then,  assuming  a  flat  input  noise  power  spectral  density,  from  Equations 
(36)  and  (38)  we  have 

Snr(«)  =  kT^F  (39) 

at  the  input.  Using  H(ju>)  as  defined  by  Equation  (32),  the  output 
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noise  power  is  given  by 


pny  =  1  Snr/|H(«)|2d«. 
2tt  -® 


(40) 


If  we  assume  toc»2^r  in  Equation  (32),  we  can  approximate  the  noise 
power  by 


ny 


5  i  • 


(41) 


By  a  simple  change  of  variables,  this  integral  can  he  put  into  the 


familiar  form  of  sinc(x),  i.e., 
P 


ny  -  I  W|£2T2j[2/T]/  s1nc2x  dx. 


(42) 


Evaluating  this  integral  with  the  restriction  of  u>c>2it,  we  have 


oo 

Pn„  -  I  Snrrc2n/sinc2x  dx  =  SnrC2T.  (43) 

ir  L  V  —  ? 

Using  Equation  (39), the  output  noise  power  can  be  obtained  as 

pny  «  kT0FC2T.  (44) 

4 

Now,  consider  an  input  signal  with  an  unknown  phase 

r(t)  =  Acos(u)ct+0)  (45) 

applied  to  the  quadrature  receiver  of  Figure  3.  It  can  be  shown  by 

reasoning  similar  to  that  which  led  to  Equation  (44)  that  the  noise 

power  at  the  output  is  given  by 

Pn  =  2FkT0C2T  , 

T 

and  that  the  desired  signal  power  at  the  output  of  the  system  Is  again 
given  by  Equation  (35).  Therefore,  we  have 
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E  =  P«  =  [ACT/2]2  =  A2T 

N  Pn  2FkTf)CzT  2kFTn 

T 

=  SBifT 
N 

E  =  [S/N]TBif  (46) 

f 

where  (S/N )  is  given  by  Equation  (13).  Equation  (46)  is  the  required 
signal-to-noise  expression  for  a  quadrature  receiver  after  signal 
processing  has  been  employed. 

4.  Pulsed  Source  with  Incoherent  Pulse  Integration. 

This  configuration  represents  a  pulsed-source  range-gated  radar. 
Such  a  radar  can  be  used  as  a  transmlssometer  by  measuring  the  reflec¬ 
tion  from  a  calibrated  target  such  as  a  corner  reflector. 

At  the  receiver,  the  signal  Is  heterodyned  to  form  an  intermediate 
frequency  signal,  passed  through  a  detector  and  sampled  sequentially  by 
range  gates.  The  purpose  of  each  range  gate  is  to  quantize  the  radar 
return  into  range  intervals.  Each  range  gate  opens  in  sequence  just 
long  enough  to  sample  the  signal  voltage  corresponding  to  a  different 
range  interval  along  the  propagation  path.  In  general,  the  range  gate 
acts  as  a  switch  (or  gate)  which  opens  and  closes  at  a  proper  instant 
of  time.  These  gates  are  activated  once  for  each  pulse-repetition  in¬ 
terval  and  produce  a  series  of  pulses  of  constant  amplitude  at  the 
gate  output  when  using  a  stationary  target. 

The  output  of  the  range  gates  is  passed  through  a  circuit  called  a 
box-car  device  which  integrates  (sums)  the  pulses.  Basically,  the  box¬ 
car  device  is  a  sample  and  hold  circuit  which  stretches  each  pulse  in 
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time  to  cover  the  entire  repetition  period.  The  purpose  of  using  a 
box-car  device  is  to  aid  in  the  filtering  and  detection  process  by 
emphasizing  the  fundamental  of  the  modulation  frequency  and  eliminating 
harmonics  of  the  pulse  repetition  frequency. 

The  receiver  system  now  examined  uses  incoherent  (post-detection) 
integration  rather  than  coherent  (pre-detection)  integration.  In  post¬ 
detection  integration,  a  smaller  signal -to-nolse  ratio  exists  when 
compared  to  pre-detection  integration.  This  is  due  primarily  to  the 
non-linear  nature  of  the  detector  which  creates  a  rectifying  effect  and 
causes  harmonic  distortion  that  decreases  the  available  signal  power. 

If  pre-detection  were  used  and  ' n 1  pulses  of  the  same  signal-to- 

noise  ratio  were  integrated,  then  the  resulting  signal -to-noise  power 

ratio  would  be  identically  'n‘  times  that  of  a  single  pulse.  However, 

for  post-uetection  integration,  a  smaller  signal -to-noise  power  ratio 

would  be  found.  In  fact,  we  can  account  for  the  reduced  signal -to-noise 

ratio  after  integration  by  using 

E  =  nEi (nifsl  (47a) 

IT  [TT 

where  E-j  (n )  is  the  integration  efficiency  factor  defined  by  Marcum  [6] 
as 

Ei (n)  =  [S/NUi  .  (47b) 

n[S/N ft 

Here,  [S/N]i  is  the  signal-to-nolse  ratio  of  a  given  pulse  required  to 
yield  a  specified  probability  of  detection,  and  [S/N]n  is  the  value  of 
signal-to-nolse  ratio  of  a  single  pulse  required  to  achieve  the  same 
probability  of  detection  when  'n'  similar  pulses  are  integrated.  E^(n) 
has  been  calculated  by  Marcum  [6]  and  plotted  by  Skolnik  [7].  The  plots 
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shown  In  Figure  5  are  for  the  integration  Improvement  factor  1=  nE-j  (n ) 
and  integration  loss  L  =  101ogio[l/Ei (n)]  both  versus  'n',  the  total 
number  of  pulses  integrated. 

For  the  present,  consider  [S/N]n  in  Equation  (47b)  to  be  equal  to 
[S/N]  in  Equation  (47a).  We  see  when  ' n '  such  pulses  are  integrated  to 
give  the  same  probability  of  detection  as  one  pulse  with  a  signal-to- 
noise  ratio  of  [S/N]j  this  value  is  approximately  the  same  as  the 
total  signal -to-noise  ratio  of  n  pulses,  [E/N]  in  Equation  (47a).  In 
general,  the  signal -to-noise  expression  with  ' n '  integrated  pulses  can 
be  given  by  -2aR 

“ITT 

[S/N]n=  P^GtGpX2cj»Lr  10  (48) 

To  use  this  equation  it  is  necessary  to  use  a  set  of  curves  as  shown  In 
Figure  5  for  each  value  of  'n'.  Therefore,  combining  Equations  (47b) 
and  (48)  yields  -2aR 

CE/N]  ■  ^  (49) 

where  Ej(n)  comes  from  Figure  5.  Note  in  this  figure  that  nf  is  the 
false  alarm  number  and  relates  to  false  alarm  probability  by  nf=l/Pfa. 
This  will  be  Important  in  probability  of  detection  calculations  to  be 
discussed  in  Section  III. 

5.  Pulsed  Source  with  £  Quadrature  Recei ver. 

This  configuration  is  a  true  pulsed  radar  which  uses  a  range  gated 
quadrature  receiver.  That  is,  the  Input  to  Figure  4a  is  given  as 

s(t)  =  Cp(t  )cosu>ct  (50) 

Here  p(t)  Is  a  unit  amplitude  pulse  train  of  width  t  and  repetition  per 
iod  Tr  contained  in  the  time  interval  (0,T).  The  output  Y(T)  of  Figure 
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4a  is  then 


T 

Y(T)  =  /  s(t)r(t)p(t)  dt  =  ACTd  (51) 

o  ~2 

where  d  Is  the  duty  cycle  of  the  signal,  described  as 

d  =  x/Tr  .  (52) 

This  integral  for  Y(T)  can  be  used  to  define  an  expression  for  the 
signal  power  by  a  method  analogous  to  Equation  (35).  That  is, 

Ps  =  [ACTd/212  .  (53) 

Once  we  choose  convenient  values  for  x,  Tr,  and  T,  the  impulse 
response  of  the  equivalent  filter  of  Figure  4b  is 

h(t)  *  Cp(t)cosu^t  (54) 

in  order  to  match  the  input  signal  s(t). 

If  the  filter  were  matched  only  to  one  pulse  of  duration  x,  we 

could  use  Equation  (44)  for  noise  power  at  the  output,  i.e., 

Px  =  kT0FC2x  .  (55) 

T 

If  we  define  the  number  of  integrated  pulses  as 


n  =  T/Tr  ,  (56) 

then  the  noise  power  of  the  filter  output  signal  with  Impulse  response 
given  by  Equation  (54)  is 


kTgFC2xn 

J 


(57) 


Therefore,  the  signal -to-noise  ratio  at  the  system  output  is 
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E  =  Pc  =  [ACTT/2Trl2  =  [A/2]2T  t  . 

*  nrr^iT7T7T  ~nToi~ 

T  7 

Since  the  power  of  a  sine  wave  of  amplitude  A  into  a  1-ohm  load  is  given 
by  S=A2/2,  then  the  processed  signal-to-noise  ratio  for  a  single 
coherently  detected  channel  is 

E  =  2STdBif 
IT  “  N 

where  [S/N]  is  defined  by  Equation  (13). 

However,  in  practice  we  will  not  know  the  phase,  so  we  must  use  a 
quadrature  receiver. 

Note  that  the  quadrature  receiver  produces  a  signal-to-noise  ratio 
at  the  output  that  is  just  one-half  as  large  as  the  phase-coherent 
detector  of  Figure  4,  this  occurs  because  the  2-phase  detectors  have 
one  channel  processing  one  component  of  signal  and  its  associated  noise 
component,  and  the  other  channel  processing  the  quadrature  signal  and 
its  associated  quadrature  noise  component.  Therefore,  we  end  up  with 
the  effects  of  the  received  signal  and  two  noise  components,  which 
halves  the  output  signal-to-noise  ratio,  i.e.. 


gives  the  signal-to-noise  ratio  at  the  system  output.  Note  this  equation 
is  very  similar  to  Equation  (46)  except  we  are  including  a  duty  cycle 
for  the  pulse. 


6.  Pulsed  Source  with  Linear  FM  Rectangular-Envelope  Pulse  Compression. 


This  transmission  configuration  is  the  pulsed  source  range-gated 


radar  discussed  in  Section  II. B. 5.  However,  the  transmitted  signal  is 
on  during  most  of  the  entire  repetition  period  but  swept  in  frequency. 

At  the  receiver,  after  the  signal  is  heterodyned  to  an  intermediate  fre¬ 
quency,  it  is  passed  through  a  matched  pulse  compression  filter  and 
detected  as  shown  in  Figure  6  [8]. 

If  the  carrier  frequency  of  a  transmitted  pulse  is  linearly  swept, 
a  puise  compression  filter  can  be  used  to  delay  one  end  of  the  pulse 
relative  to  the  other.  At  the  filter  output,  this  produces  a  narrowed 
pulse  of  greater  peak  amplitude.  The  linear  time-delay  characteristic 
of  the  filter  acts  to  delay  the  high-frequency  components  at  the  start 
of  the  input  pulse  more  than  the  low-frequency  components  at  the  end  of 
the  pulse.  The  frequency  components  experience  a  porportional  delay  so 
that  the  net  result  is  a  time  compression  of  the  pulse.  The  effects  of 
pulse  compression  are  summarized  in  Figure  7. 

It  can  be  shown,  as  in  Appendix  8,  the  pre-detection  signal-to- 
noise  ratio  for  the  radar  transmissometer  is 


S  =  PtKGtGrltLroX2  10 
N  ' 


-2aR 

~1U 


(59) 


where  K  is  the  pulse  compression  ratio  defined  by 


K  =  T/t  . 


(60) 


where  T  and  x  are  defined  in  Figure  7.  The  ratio  of  pulse  widths 
before  and  after  pulse  compression  is  shown  in  Figure  o. 

If  we  consider  a  pulse  compression  radar  with  a  quadrature  recei¬ 
ver  system  as  in  Section  II.  R.5,  then  the  output  signal -to-noise  ratio 
can  be  found  from  Equation  (58)  to  he 
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Figure  6.  Typical  Intrapulse  Linear  Frequency  Pulse  Compression  Radar 
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Figure  7.  Idealized  Pulse  Compression  Characteristics. 

(a.)  Wide-pulse  envelope,  (b.)  Carrier  freq¬ 
uency  modulation,  (c.)  Filter  time-delay 
characteristics,  (d.)  Compressed  pulse  enve¬ 
lope,  (e.)  Input-output  waveforms  of  compres¬ 
sion  filter. 
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-2aR 

TT 

E  =  PtKGtGrLtLraX^TdRif  10  .  (61) 

F  t(4.)¥k^ — 


D.  Summary 

In  this  section,  signal -to-noise  ratio  expressions  before  signal 
processing  were  derived  in  Section  II, A  for  both  double-ended  and 
single-ended  systems,  they  are  given  by  Equations  (7)  and  (13). 

Also,  post-detection  signal -to-noise  ratio  expressions  for  the  si 
systems  were  derived  in  Section  II.B.  These  signal -to-noise  ratio 
expressions  are  given  by  Equations  (16),  (46),  (49),  (58),  and  (61). 

In  all  cases,  the  pre-detection  signal -to-noise  ratios  to  be  used  are 
tnose  derived  in  Section  II. A. 

Pre-detection  signal -to-noise  ratios  are  denoted  by  S/N  and  post¬ 
detection  signal -to-noise  ratios  denoted  by  E/N. 

Note  that  all  the  equations  derived  in  this  section  are  dependent 
on  the  attenuation  parameter  a.  A  derivation  of  equations  for  a  will 
be  given  in  Section  IV  for  attenuation  in  rain  and  fog. 
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SECTION  III 


DERIVATION  OF  PROBABILITY  OF  DETECTION 
For  transmissometers,  the  output  (post-detection)  signal -to-noise 
ratio  is  probably  the  most  useful  specification  of  performance.  For 
radars,  most  radar  designers  prefer  to  use  a  specification  of  probabil¬ 
ity  of  detection  and  probability  of  false  alarm  rather  than  signal-to- 
noise  ratio.  Thus,  in  this  section  a  derivation  of  expressions  for  the 
probability  of  detection  as  a  function  of  signal -to-noise  ratio  will  be 
presented.  Since  the  signal -to-noise  ratio  is  directly  related  to 
range  as  shown  by  the  equations  in  Section  II,  it  will  be  possible  to 
derive  probability  of  detection  expressions  in  terms  of  range  for  given 
system  parameters  and  meteorological  conditions.  In  addition,  a  speci¬ 
fication  of  false  alarm  probability  must  be  given.  The  discussion  pre¬ 
sented  in  this  section  is  related  closely  to  the  modified  Neyman -Pearson 
detection  theory  described  by  Skolnik  [91. 

The  method  presented  in  this  section  is  derived  with  the  assumption 
of  either  a  square-law  or  linear  detector  since  the  effect  on  the  de¬ 
tection  probability  by  assuming  one  instead  of  the  other  is  small.  We 
will  also  apply  this  technique  to  quadrature  detection  systems  since  the 
probability  of  detection  will  be  only  a  function  of  processed  signal -to- 
noise  ratio  and  false  alarm  probability. 

A.  Description  of  Noise  in  the  Radar  System. 

To  successfully  define  the  probability  of  detection  in  terms  of 
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signal -to-noise  ratio  (or  range)  we  need  a  suitable  representation  of 
the  noise.  One  of  the  best  noise  models  is  discussed  by  Skolnik  where 
we  consider  narrowband  Gaussian  noise  which  is  both  amplitude  and  phase 
modulated.  For  instance,  if  we  let  n(t)  represent  the  random  noise 

IN/ 

function  (where  ~  denotes  a  random  variable),  then  set 

n(t)  =  g(t)cos[>ct+4>(t)]  (62) 


where  g(t)  and  <j>(t)  are  random  variables  representing  amplitude  and 
phase  which  vary  slowly  when  compared  to  the  cosine  factor.  Thus, 
n(t)  =  g(t)cos<J>(t)cosioct  -  g(t)sin<j>(t)sinu)Ct  .  (63) 


However,  these  leading  terms  are  simply  time  varying 
To  simplify  the  notation,  consider 

a(t)  =  g(t)cos<Kt) 

<w 

and 

h(t)  =  g{t)sin<j>(t). 


random  variables. 

(64) 

(65) 


so  that 

1/2 

g(t)  =  [a2(t)+b2(t)]  (66) 


.  b(t) 

4>(t )  =  tan-1  ~  .  (67) 

"  a'(t) 


Therefore,  the  noise  function  can  be  represented  as 

n(t)  =  a(t)cosoJCt-b(t)sinu)Ct  .  (68) 


From  a  careful  examination  of  the  statistical  properties  of  typical 
receiver  noise,  it  has  been  found  that 

(1.)  The  random  variables  a(t)  and  b(t)  are  statistically 
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independent.  Since  they  are  uncorrelated,  it  follows  that 
E[a(t)b(t)]=0. 

(2.)  The  mean  values  E[a(t)]=E[b(t)]  are  zero,  and  the  variances 
E[a^(t)]=E[b2(t)]  are  equal  to  a  constant  <|>o. 

Therefore,  we  can  write  the  signal  input  to  the  detector  in  the  form 
y(t)  »  [f  (t)+a(t)]cosa)ct-[b(t)]sina)Ct  (69a) 

r+s  m 

if  a  signal  is  present,  and 

y(t)  =  a(t)cosojct-b(t)sina)ct  (69b) 

rv  /V  (V 

if  no  signal  is  present.  For  a  pulsed  signal, 

A  |t|<  T/2 

f(t)  ={  (70) 

0  |t|>  T/2 

represents  the  constant  amplitude  of  a  pulse.  For  a  CW  signal,  A///T 
represents  the  r.m.s.  amplitude  of  the  carrier  when  f(t)=A  is  used. 

If  y(c)  is  passed  through  a  noncoherent  detector  (such  as  a  half  or 
full-wave  rectifier  followed  by  a  low  pass  filter,)  then  all  phase 
information  is  lost  in  the  detection  process.  The  output  of  the 
detector,  r(t),  assuming  unity  gain  in  the  detection  process  will  be 
given  by 

r(t)  =  ([f(t)+a(t)]2+[b(t)]2}l/2  (71) 

if  a  signal  is  present,  and 

r(t)  -  {ra(t)]2*[>(t)32)l/2  (72) 

/N*  TV 

if  no  signal  is  present. 


B.  Derivation  o£  the  Probability  of  Detection. 

For  our  work,  we  are  primarily  Interested  in  Type  II  error.  This 
error  occurs  during  the  detection  process  if  we  decide  no  signal  is 


present  when  it  actually  is  present.  If  we  choose  an  arbitrary 
threshold  level  Vf,  the  probability  of  having  a  Type  II  error  can  be 
defined  as 


VT 

Pod  =  I  P(r(t))  dr  (73) 

0 


and  the  probability  of  detection  defined  as 


oo 

Pd  -  /  P(r(t))  dr  (74) 

VT 

where  p(r)  is  the  probability  density  function  of  r(t),  which  will  be 
assumed  to  be  a  Gaussian  random  variable.  Me  can  describe  the  probabil¬ 
ity  of  detection  and  probability  of  Type  II  error  as  shown  in  Figure  9. 
Also,  for  convenience  in  the  notation,  let's  consider 

z(t)  =  f  (t)+a(t)  (75) 


in  Equation  (69a).  Since  the  random  variables  a(t)  and  b(t)  are  statis- 
tically  independent,  it  follows  that  z(t)  and  b(t)  are  also  independent. 

M  (V 

Therefore,  Equation  (69a)  becomes 

r(t)  -  l[z(t)]2+[b(t)]2(l/2  (76) 

iv 

when  a  signal  is  present. 

Now,  consider  the  case  where  f(t)=A  for  a  pulsed  radar.  Then  z(t)= 
A+a(t)  is  a  Gaussian  random  variable  with  mean  A.  That  Is, 

TV 

E[z(t)]  =  E[A]  +  B[a(t)J  =  A  (77) 

TV  TV 

since  E[a(t)]=0  and  the  expected  value  of  a  constant  is  the  value  of  the 
constant.  Also,  we  can  show 

E[(z(t)-A)23=E[a2(t)>i|70  .  (78) 

From  the  bottom  of  page  30,  we  see  E[b2(t)]]=  \|>o.  Since  we  have  shown 
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E[z(t)]=  A  and  E[b(t)]=  0,  we  can  write  the  probability  density  func¬ 


tions  for  z (t )  and  b(t)  as 


p[z(t)]  =  1  e 


tz{t)-A]2' 


b2(t) 


p[b(t)]  =  1  e  . 

~  TTrnj^ 

However,  since  a{t)  and  b(t)  are  uncorrelated,  we  can  show  z(t)  and  b(t) 

<v  /V  ^  ^ 

are  independent,  i.e., 

E[z(t)b(t)]  =  E[(a(t)+A)b(t)]  =  E[a(t)b(t)]  +  E[Ab(t)1 


=  AE[b(t)] 


E[z(t)b(t)]  =  0  . 

/s/  rs* 

Knowing  z(t)  and  b(t)  are  independent,  we  see  the  joint  probability 
density  function  can  be  expressed  as  the  product 

p[z(t)b(t)]  =  p[z (t)]p[b(t)]  . 

/V  ^ 

That  is,  the  joint  probability  density  function  becomes 

frz(t)-A]2+b2(t)l 


p[z(t)b(t)]  =  J_  e  .  (79) 

~  ~  2inJ>o 

We  wish  to  convert  this  expression  Into  a  density  function  in  terms  of 
r(t)  at  the  output  of  the  noncoherent  detector. 
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noise  only 


:  Pd=  Probability  of  Detection 


:  Pfa*  Probability  of  Falsa  Alarm 


Figure  9.  Probability  of  Detection  and  Probability 

of  False  Alarm  In  terms  of  density  function 
plot  p(r). 


Suppose  we  consider  the  two  random  variables  z(t)  and  b(t)  as 

/V 

rectangular  coordinates.  Then  we  are  able  to  convert  to  polar  coordi¬ 
nates  where  computations  are  much  easier.  This  is  equivalent  to  using 
the  polar  representation  of  the  noise  given  by  Equation  (66)  and  (67). 
By  Equation  (79),  we  actually  mean 

-  [f- A]2+62] 

-  - 

p[z-dz<z(t)<z+dz  ,  6-dfj<b(t)<6+d6l  =  1  e  dz  dti.  (80) 

L  T  -  T  T  ~  "7J 


This  defines  the  probability  of  finding  z(t),  b(t)  in  the  small  neigh- 

~  /V 

borhood  z,b  in  a  cartesian  system.  We  can  transform  to  polar  coordi¬ 


nates  by  using 


z  =  rcosQ  b  =  rsinQ 


so  the  elemental  area  dzdb  becomes  rdrde.  These  coordinate  transforma 
tions  are  sketched  in  the  drawing  below. 


de\  , 


Therefore,  we  can  write  , 

-  [rcos 6-A]2+[rsine]2 
_  _  _  _  _  _  _  Z*0  J 

p[>-dr<r<r+dr  ,  8-d9<9<9+d0l  =  1  e  rdrd9 


or 


p(r,9)  =  r  e 


=  r _ e 

2**o 


-  f[rcos9-A]2+[rsin9]2 
L  2*o  ’ 


r2+A2_2rAcos9 

?*5  ; 


(81) 


Now,  the  total  probability  of  finding  r  between  r+dr  and  r-dr  corre- 

2  T 

sponds  to  a  ring  of  width  r  about  the  circle  r=r.  Therefore, 


2*  2* 
p(r)  =  /  p(r,9)d©  «  /  r  e 
0  ~  ~  0  2**o 


2* 


-  jr?+»*  -2rAcos9 


2*0 


Jd9. 


r^+A2 

A r cos 9j 

2*0  . 

2* 

«  r 

*0  J 

d9.  (82) 


0  -  *0  2*  0 

This  latter  integral  is  defined  as  the  modified  Bessel  function  of  the 

first  kind  and  order  zero  with  argument  Ar/*o.  That  is, 

r2+A2] 

’Tio" 


p(r)  =  r  e 
*0 


Jo[Ar/*ol  . 


(83) 


Now,  knowing  the  probability  density  function  p(r),  we  can  calcu¬ 
late  the  probability  of  detection  via  Equation  (74). 

fr2+A2' 


pd  =  /p(r)dr  =  f  re 
Vy  Vj  iJ>q 


rwj 


JoCAr/*o]dr  .  (84) 


This  integral  is  evaluated  by  Skolnlk  [9]  in  terms  of  an  Infinite 
sum,  so  that 
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F-  ■■■  ■  jmmim 
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kv 


(Vt-A)2 


Pd  =  rer«] + 


♦o 


2A  |2tt 

l/*0 


(Vt-A)2 
Vt-a+  1+  4»n 
"W - 5^ - 


<Po 


(85) 


Note:  erf(z)  is  the  error  function  defined  as 

z  -u2 

erf  (z)  =  _2_/e  du 
/¥  0 


(86) 


A  valid  series  approximation  to  Equation  (85)  can  be  found  when 
we  choose  a  sensible  threshold,  i.e.,  (Vj-A)2/i|>o<1  and  we  have  a  good 
signal-to-noise  ratio  so  A2/ij>g»l  .  If  these  conditons  are  valid,  as 
often  occurs  for  fairly  large  signal-to-noise  ratios,  we  can  use 


Pd 


r 


-er  f/V-r-AM  +  e 


_  (Vt-a)2 


2A/f? 
V  *0 


(Vr-A) 2 
1-  Vt-A  +  1+  M _ 


~w 

<J>0 


(87) 


According  to  Skolnik,  if  we  examine  the  envelope  of  the  IF  sine 
wave  rather  than  the  signal  after  detection,  it  is  then  possible  for 
the  equation  for  probability  of  detection  to  be  converted  to  a  power 
relationship  by  defining  signal-to-noise  ratio  at  the  output  of  the 
detector  as  E/N.  That  is, 

r 

[A//^1]2  for  a  pulsed  signal 
[A/ZZlj/ol2  for  a  CW  signal. 

Using  these  relationships,  we  can  see 

-UM- 


E  =  Power  in  Signal  = 
TT  Power  in  Noise 


(88) 


Pd  =  lfl-erfG^  )~1+ 


Am)  j  /  8tt[E/N]'  v 

m 

where  [E/N]>> (A/Vj)2  for  a  pulsed  signal  and  [E/N]»(A/2Vy)2  for  a  CW 
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.V. 


signal.  We  will  see  in  Section  V  that  Equation  (89)  proves  to  he  accur¬ 
ate  for  signal -to-noise  ratios  in  excess  of  0  dB. 

C,  Derivation  of  Probability  of  False  Alarm. 

The  probability  of  false  alarm  is  the  probability  that  the 
threshold  voltage  is  exceeded  when  no  signal  is  present.  Commonly 
called  Type  I  error,  the  false  alarm  probability  is  defined  by  the 
integral 

00 

pfa  =  Jp(r)dr  (90) 

VT 

where  p(r)  is  given  by  Equation  (83)  when  no  signal  is  present.  Because 
t^Ar/<|>o)=l  when  there  is  no  signal,  we  find 

-p2 

2+0 

p(r)  =  r  e  (91) 

*0 

This  probability  density  function  is  commonly  known  as  the  Rayleigh 
density  function. 

Substituting  Equation  (91)  into  (90)  and  integrating  gives 

-pi 

Pfa  »  e  (92) 

for  false  alarm  probability.  Using  Equation  (89)  allows  a  definition 
of  probability  of  false  alarm  in  terms  of  signal -to-noise  ratios  at  the 


output  of  the  detector.  That  is, 


I  Vt2[E/N] 
1  "  ft?- 


pfa 


I 


for  pulse  signal 


Vtz[E/N] 
-  J  -jr 


(93) 


for  CW  signal 


Solving  this  equation  for  the  ratio  of  V t/A,  yields 


(94) 


I 


A 


I 


r. 

i 

i 
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for  pulse  signal 
for  CW  signal. 


Substitution  of  Equation  (94)  into  Equation  (89)  allows  us  to  obtain 
a  probability  of  detection  in  terms  of  a  specified  false  alarm 
probability  and  processed  signal -to-noise  ratio. 

The  processed  signal -to-noise  expressions  (E/N)  have  been  derived 
in  Section  II. B.  as  a  function  of  system  parameters.  Therefore,  once 
the  radar  system  designer  has  specified  a  probability  of  false  alarm, 
the  required  threshold  ratio  can  be  determined  from  equation  (94)  and 
the  probability  of  detection  from  Equation  (89).  Or,  as  often  is 
needed,  the  radar  designer  may  desire  to  specify  both  probability  of 
detection  and  a  probability  of  false  alarm.  In  this  case,  the  required 
signal -to-noise  ratio  at  the  detector  output  can  be  calculated  and  the 
range  at  which  it  is  equalled  or  exceeded  can  be  found.  Both  methods 
will  be  discussed  with  graphs  and  computer  algorithms  in  Section  V. 

Typical  plots  showing  probability  of  detection  as  a  function  of 
signal -to-noise  ratio  or  range  are  shown  below  in  Figures  10  and  11.  In 
both  cases  a  specification  of  false  alarm  probability  must  be  given. 
Figure  10  can  be  plotted  by  using  Equation  (89).  Figure  11  is  plotted 
by  using  Equation  (89)  and  the  signal -to-noise  ratio  expressions  as  a 
function  of  range  derived  in  Section  II.  Therefore,  for  all  of  the 
radar  systems,  we  can  plot  probability  of  detection  either  in  terms  of 
signal -to-noise  ratio  or  range  for  a  fixed  probability  of  false  alarm. 

Finally,  the  proper  choice  of  the  signal -to-noise  ratio  given  by 


•l- 


JL>  JLm  fc 
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Equation  (88)  will  be  given  for  each  system  discussed  in  Section  II. 

The  derivation  of  the  equations  for  probability  of  detection  deals 
only  with  the  envelope  of  the  IF  output  signal,  which  is  a  sine  wave. 
Therefore,  A  is  the  amplitude  of  the  IF  signal  in  the  absence  of  noise 
and  S=a2/2  is  the  power  in  the  IF  channel  due  to  signal  when  a  pulse  is 
present. 

For  the  two-way  transmission  and  square-law  detection  system  of 
Section  II. B. 2,  we  use  a  "chopped"  signal  with  square-law  detector 
followed  by  an  audio  filter.  In  Equation  (16),  we  have  S  representing 
the  carrier  power  at  the  output  of  the  intermediate  frequency  amplifier 
and  E  representing  Ep2/2  where  Ep  is  the  peak  value  of  the  audio  sine 
wave.  If  we  consider  that  this  audio  wave  will  be  envelope  detected, 
then  this  reflects  Equation  (88)  be  used  for  a  CW  signal  and  E=A2/2. 
Therefore,  in  S/N  for  this  system,  interpret  P^  as  the  carrier  power  so 
for  "chopped"  signal  and  square-law  detection  use  the  bottom  right  side 
of  Equations  (88)  and  (94). 

For  the  CW  signal  transmission  with  a  quadrature  receiver  of  Sec¬ 
tion  II. 8. 3,  the  post-detection  "power"  was  obtained  by  taking  a  sample 
of  signals  with  mean  value  of  A1?,  i.e.,  E=a2.  Therefore,  in  S/N  for 
this  system,  interpret  P^  as  the  carrier  power  which  is  equal  to  the 
average  or  peak  power.  In  this  case,  use  the  top  right  side  of  Equation 
(88)  to  determine  post-detection  signal -to-nolse  ratio  E/N. 

For  the  pulsed  source  system  with  Incoherent  pulse  Integration 
discussed  in  Section  II. B. 4,  the  pre-detection  S/N  refers  to  a  sine  wave 
signal  and  P^  represents  the  peak  power  radiated  during  a  pulse.  If  A 
is  the  envelope  of  the  IF  signal  in  the  absence  of  noise,  then  S=A^/2 
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is  the  power  in  the  IF  channel  due  to  signal  when  a  pulse  is  present. 

In  this  case,  use  the  bottom  right  side  of  Equations  (88)  and  (94). 

For  the  pulsed  source  with  quadrature  receiver  of  Section  II. 8.5 
and  the  pulsed  source  with  linear  FM  rectangular  envelope  pulse  compres¬ 
sion  of  Section  II. 8. 6,  the  same  reasoning  applies  to  E/N  as  to  the  CW 
signal  and  quadrature  receiver  discussed  above.  In  each  of  these  cases, 
the  post-detection  "power"  was  taken  from  a  sample  of  signals.  That  is, 
Pt  is  the  peak  power  transmitted  during  a  pulse  so  that  S=A^/2  (i.e., 
S=(a2cos0)/2  in  one  channel  and  S=(A2sin6)/2  in  the  other  channel)  leads 
to  Equation  (58).  In  the  pulse  compression  case,  the  peak  and  average 
power  are  the  same  if  pulse  width  before  pulse  compression  T,  equals  the 
repetition  period  Tr.  Once  again,  the  peak  power  transmitted  during  a 
pulse  is  S=A2/2.  For  these  two  cases,  we  use  the  top  right  side  of 
Equations  (88)  and  (94)  for  the  probability  of  detection  calculation. 
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SECTION  IV 


ATTENUATION  DUE  TO  RAIN  AND  FOG 

The  attenuation  parameter  o  in  the  signal -to-nolse  ratio  expres¬ 
sions  is  highly  dependent  upon  weather  conditions  and  frequency.  Ap¬ 
proximate  values  at  sea  level  in  clear  air  are  shown  In  Figure  1,  but 
they  are  subject  to  variation  because  the  absorption  by  oxygen  and  water 
vapor  changes  substantially  with  weather  conditions.  Even  more  pro¬ 
nounced  is  the  attenuation  due  to  rain  and  fog  for  which  we  now  wish  to 
present  a  model . 

For  our  path  attenuation  model,  assume  that  rain  and  fog  consist 
of  small  spherical  water  droplets  of  complex  refractive  Index  to  account 
for  dissipation  and  scattering  effects.  The  scattering  effects  from  a 
sphere  were  examined  in  detail  by  Mie  [10], are  presented  in  an  excellent 
tutorial  book  by  Kerker  [11],  and  have  been  recently  utilized  in  propa¬ 
gation  algorithms  by  Ruyukdura  [2].  Therefore,  a  detailed  examination 
of  these  scattering  effects  will  not  be  given  in  this  thesis.  Instead 
the  work  by  Buyukdura  will  be  modified  to  include  effects  of  frequency 
and  temperature  on  the  complex  index  of  refraction  for  water.  These 
effects  are  discussed  in  an  article  by  Peter  S.  Ray  [12].  Corrections 
to  this  article  are  given  by  [13]  and  the  calculations  for  the  complex 
Index  of  refraction  In  terms  of  frequency  and  temperature  are  shown  In 
Appendix  A. 
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In  general,  the  attenuation  factor  caused  by  a  spherical  drop-size 
distribution  using  the  extinction  efficiency  factor  is  given  by  Buyuk- 
dura  [2]  as 

00 

a  =  4.343  /  QEXT(a)N(a)da  (95) 

0 

where  a  is  specific  attenuation  in  dB/km,  a  is  droplet  radius  in  meters, 
N(a)  is  the  volume  density  of  droplets  with  radius  a  in  units  of 
and  QgxT  is  the  extinction  efficiency  factor 

oo 

Qfxt  -  -2  I(2n+l)[Re(an+bn)].  (96) 

T&OFP  n=l 

Here  kg  is  the  wave  number  in  free  space,  and  an,  bn  are  the  spherical 
Mie  coefficients  defined  in  Equations  (A63)  and  (A64)  of  Buyukdura  [2]. 
A,  Oropsize  Distribution  for  Rain. 

For  our  work,  we  will  use  a  Marshal  1 -Palmer  dropsize  distribution 
of  the  form 

N(a)  =  N0e-2Aa  (97) 

where  A  =  a^R^mm-1), 

Nq  =  8000  [drops/(m3«mm)], 

®1  =  4.1, 

6  =  0.21, 

and  R  is  the  rain  rate  In  mm/hr,  a  Is  the  droplet  radius  in  mm.  There¬ 
fore,  substituting  Equations  (96)  and  (97)  into  Equation  (95)  gives 
an  expression  for  the  extinction 

00 

a  =  ag  +  4. 343/Q£XT(a )N(a )da  (dB/km)  (98) 

0 


where  the  constant  q)=0.02  dB/km  at  35  GHz,  ao=0.3  dB/km  at  94  GHz,  ag= 
1.4  dB/km  at  140  GHz,  and  oq=4  dB/km  at  240  GHz  in  clear  atmosphere  [14] 
B.  Dropsize  Distribution  for  Fog. 


Dropsize  distributions  for  advection  fogs  near  seacoasts  have  been 


extensively  studied  in  the  works  of  Kunkel  [15]  and  Mallow  [16].  We 
have  found  no  detailed  analyses  for  radiative  fogs  common  in  the  midwest 


United  States.  Therefore,  the  description  of  a  fog  presented  in  this 
section  is  based  upon  a  sparse  advectlon  fog  model  derived  by  Kunkel. 

A  widely  used  dropsize  distribution  is  based  upon  Kunkel  [15] 

,W 

N (A )  =  CAa  e  y  (99) 

where  N(A)  is  the  number  density  of  droplets  of  relative  radius  A=a/Rm. 
Here,  a  Is  the  actual  radius  of  the  fog  droplet  and  Rm  is  the  drop 
radius  with  the  maximum  number  density  in  the  distribution.  The  normal¬ 


ization  constant  is  generally  determined  from  the  Integral 

oo 

N  =  Rm  /  N(A)  dA 
0  y 

oo  -g'  A 

=  RmC  /  A“’e  y  dA  (100) 

0 


which  can  be  solved  for  C  as 

.  m 

C  =  +Ny  [a  /y] 

^m 


1 

T  [a*  +1 


(101) 


where  N  is  the  total  number  of  droplets  in  the  distribution  scaled  to 


conform  to  the  total  liquid  water  content  of  a  given  fog.  Note  that  the 
(+)  sign  Is  used  when  a'and  y  are  positive.  Likewise,  the  (-)  sign  is 
used  when  a'and  y  are  negative. 

As  Equations  (99)  and  (101)  are  given,  we  must  determine  the 
scaling  factor  used  to  conform  to  the  total  liquid  water  content  of  the 
fog.  The  liquid  water  content,  0,  i.e.,  mass  per  unit  volume,  is  given 
by 
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D  =  pj  4ira3N(a)  da 
0  ~1 

where  a  Is  the  actual  droplet  radius  and  p  is  the  density  of  water,  then 


D  =  pRm4/  4ttA3N(A)  dA 
0  T 


(102) 


we  see 


R ml  N(A)  dA 

N  =  0 _ 

IT  « 

pRm4/  4irA3N(A)dA 
0  T 


(103) 


0Rm/  N(A)  dA 

N  =  _ 0 _ 

00 

pRm4/  4irA3N(A)dA 
0  T 


(104) 


Substituting  Equation  (99)  and  integrating  gives  N  as  a  function  of 

liquid  water  content  and  relative  droplet  radius. 

«b  -a  A 
DRm/CA“  e  y  dA 

N  =  0 _  (105) 

m  -a1 AV 

4irpRm4/  A(3+a')e*T“dA 

T  0 


which  can  be  solved  in  closed  form  as 


(3/y) 

N  =  DRm  [o'/y]  r[(l+a'  )/y1  • 

4*pRnT  r[(4+et'  )/y] 


Substituting  Into  Equation  (101)  and  then  into  Equation  (99)  shows 

,  ,  C(4+«*)/y]  ,-a  AY 

N(A,0)  =  T  D|yI  [a'/y] _  A®  e  y  (106) 

TirpR^4  rt(4+a'/y] 

7 
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for  the  drop-size  distribution  in  fog.  Mallow  has  examined  several  fogs 

and  has  determined  a',  y,  and  Rm  for  a  light  advection  fog  to  be 

Rm  =  15  ym 
a1  =  -2.20 
y  =  -4.54 

which  yields  units  of  N(A,D)  in  droplets/(ym)4. 

We  intend  to  use  the  method  presented  by  Buyukdura  to  determine 
the  attenuation  factor.  This  requires  the  units  of  N(A,D)  to  be  in 
droplets/m3*mm  for  D  and  p  in  g/m3  and  a  in  ym.  Therefore,  we  see 

,  ,  C(4+a)/y]  -aAY 

N(A,0)  =  T  d!yI1021  [a/y]  A«e  y  (107) 

TW - r[(4+a)/y]' 

1 

in  drop1ets/(m3*mm).  Substituting  a',y,  and  Rm  with  the  knowledge  that 
A=a/Rm  and  p=106  g/m3  gives. 

5  -4.54 

N(a,n)  «  2.9536602xl0l°Da-2*2e-1»0588176xl0  a  ’  (108  ) 

in  droplets/(m3*  mm)  for  a  in  ym. 

We  will  allow  D  to  vary  between  0.021  and  0.148  g/m3  and  a  to  vary 
between  10  and  60  ym  for  a  useful  algorithm  to  calculate  attenuation  by 
fog.  Equation  (108)  will  be  substituted  Into  Equation  (98)  for  the 
calculation  of  the  extinction  factor  for  fog.  Finally,  this  model 
does  not  include  effects  of  multiple  scattering  from  each  fog  droplet. 
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SECTION  V 


RESULTS  ANO  DISCUSSION 

The  purpose  of  this  thesis  is  to  compare  the  performance  of  three 
types  of  transmissometers  and  three  types  of  radar  systems  as  a  function 
of  meteorological  conditions  in  the  millimeter-wave  windows  at  35,  94, 
140,  and  240  GHz.  For  the  radars,  a  reasonable  criterion  Is  the  range 
at  which  they  would  be  effective,  i.e.,  at  which  a  specified  high  detec¬ 
tion  probability  would  be  achieved,  subject  to  a  reasonable  false  alarm 
probability.  For  the  transmissometers,  the  corresponding  criterion  is 
the  range  at  which  a  reasonable  signal -to-noise  ratio  will  be  achieved. 
These  range  criteria  have  the  desirable  feature  that  they  do  not  depend 
strongly  on  the  detailed  system  parameters,  especially  when  attenuation 
due  to  rain  or  fog  Is  strong.  The  plots  of  range  versus  rain  rate  or 
fog  density  will  be  found  In  Figures  52  through  58;  they  are  discussed 
in  more  detail  In  Sections  V.B  and  V.C,  respectively.  Other  useful 
criteria,  such  as  signal-to-noise  ratio  and  detection  probability  as 
functions  of  range,  are  calculated  in  the  process;  these  will  be  discus¬ 
sed  and  sample  results  will  be  presented  as  they  are  encountered  in  the 
discussion. 

For  each  of  the  six  types  of  Instrumentation  discussed  in  Section 
II,  the  signal-to-noise  ratio  before  signal  processing  is  computed  as  a 
function  of  range  by  use  of  Equation  (7)  for  double-ended  systems  and 
Equation  (13)  for  single-ended  systems.  The  range  R  is  varied  from  500 
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meters  to  8  kilometers  by  increments  of  100  meters.  Then,  one  of  Equa¬ 
tions  (16),  (46),  (49),  (58),  and  (61),  as  appropriate  for  the  particu¬ 
lar  system,  is  used  to  calculate  the  signal-to-noise  ratio  after  detec¬ 
tion. 

In  our  system  calculations,  the  propagation  path  is  assumed  to  be 

uniform,  or  to  be  comprised  of  layers  each  with  homogeneous  specific 

attenuation.  The  overall  effect  of  these  layers  is  seen  by  substitution 

of  Io|Ri  for  cR  in  Equations  (7)  and  (13).  For  the  summation,  is 
i 

the  specific  attenuation  and  is  the  one-way  path  length  through  the 
the  ith  layer. 

Computer  algorithms  were  also  developed  to  calculate  and  plot  the 
probability  of  detection  as  a  function  of  processed  signal-to-noise 
ratio,  and  and  as  a  function  of  range  for  each  of  the  radar  systems. 

In  these  radar  detection  calculations,  a  specification  of  probability 
of  false  alarm  and  radar  cross  section  area  must  be  given.  These 
algorithms  use  Equations  (7)  and  (13)  to  calculate  signal-to-noise 
ratio  before  signal  processing.  Then  one  of  Equations  (16),  (46),  (49), 
(58),  or  (61)  depending  upon  the  specific  type  of  system  Is  used  to 
determine  the  signal-to-noise  ratio  after  detection  and  signal  proces¬ 
sing.  Using  the  results  presented  in  Section  III,  the  false  alarm 
probability  and  the  processed  signal-to-noise  expressions  can  be  used 
in  Equation  (94)  to  yield  Vj/A,  and  the  probability  of  detection  is 
then  computed  according  to  Equation  (89). 

Also,  computer  routines  have  been  Implemented  to  calculate  and  plot 
the  maximum  attainable  range  as  a  function  of  specific  attenuation  a, 
rain  rate,  or  liquid  water  content  In  fog.  From  the  specification  of  a 


fixed  probability  of  detection  and  probability  of  false  alarm,  Equa- 
ions  (89)  and  (94)  are  solved  for  the  processed  signal -to-noise  ratio 
by  a  search  technique  which  increases  E/N  from  a  value  of  1.0  in  steps 
of  0.001  until  the  desired  value  of  P<j  is  reached.  Then  specific  atten¬ 
uation  is  varied  from  30  to  0  dB/km  in  increments  of  -0.1  dB/km.  From 
the  signal-to-noise  ratio  expressions  before  signal  processing,  an 
iterative  solution  for  maximum  range  can  be  determined  for  large  values 
of  a  by  using  Equations  (7)  and  (13).  From  Equation  (7),  an  iterative 
solution  for  range  is  given  by 


Rn+1  -  JL2 

a 


login 


Pf-GfGr  ^Lf-Lp 
(  4  TT  )  ^  (cTqBF  (  S /N  ) 


-  ZlogioRn 


(109) 


for  n  ar.  integer.  In  this  equation,  S/N  is  the  signal-to-noise  ratio  at 
the  output  of  the  intermediate  frequency  amplifier  corresponding  to  the 
processed  signal-to-noise  ratio  for  a  fixed  probability  of  detection  and 
probability  of  false  alarm.  By  a  similar  approach,  propagation  systems 
over  a  two-way  path  can  use  Equation  (13)  to  solve  for  range  so  that 


Rn+1  =  10 
7oi 


login 


PfGfGrX^LtLrO 


-  41ogmRn 


(110) 
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In  these  iterative  solutions,  for  a=30  dB/km,  the  iterations  begin 
with  Ro=0  m.  For  all  other  values  of  specific  attenuation,  the  solution 
corresponding  to  the  preceeding  a  value  R(a+0.1)  is  used  for  R().  In  all 
cases,  the  iterations  were  terminated  when  |  Rn+i-Rn  |  <  1  m.  As  previ¬ 
ously  stated,  this  method  was  found  by  Buyukdura  to  converge  only  for 
large  values  of  specific  attenuation. 

A  different  iterative  scheme  is  used  for  small  values  of  a  when 
Buyukdura's  method  described  by  Equations  (109)  and  (110)  fails.  For 
this  case,  the  range  R  appearing  in  the  denominator  of  Equation  (7)  is 
used  on  the  left  hand  side  of  the  equation.  That  is,  for  the  two-ended 
system 


R"H  ftffiOrl 


-aRn 

10  Xi. 


(Ill) 


Similarly,  for  systems  of  the  radar  type.  Equation  (13.)  leads  to 

, - , - , 

R„+l  =  /PtGt6rX2LtLro  10  ^iT 

/  (fe?CToBr(S/iiT  .  (u 


1  (»'.TskT0HF(5/H)  .  (112) 

Again,  R(a+0.1)  was  taken  for  the  first  value  of  R  to  be  used  on  the 
right-hand  side  of  the  equation  and  the  iterations  were  terminated  when 
I  Rn+l_Rnl < 

For  all  calculations  and  graphs  presented  in  this  thesis,  the 


values  used  as  input  parameters  for  computer  algorithms  are  shown  in 


Table  I,  unless  specified  otherwise.  These  values  are  typical  of 
experimental  millimeter  wave  transmission  systems.  Note  that  all 
values  in  the  table  are  not  used  in  each  calculation.  For  some  graphs 


of  the  six  types  of  instrumentation  presented  in  this  section,  it  is 
assumed  the  propagation  path  is  through  a  homogeneous  medium  with 
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TABLE  I 


Input  Data  for  Computer  Algorithms 


Input  Parameter  Frequency  (GHz) 


35 

94 

140 

240 

Transmitted  Power*(Pt): 

200 

200 

50 

20 

mW 

Transmitting  Antenna  Gain  (G^): 

47 

47 

47 

47 

dB 

Receiving  Antenna  Gain  (Gr): 

47 

47 

47 

47 

dB 

Total  System  Loss  (L^+Lp): 

6 

8 

9 

11 

dB 

Noise  Figure  (F): 

3 

4 

7 

15 

dB 

IF  Amplifier  Bandwidth  (B-j f ) : 

10 

10 

10 

10 

MHz 

Audio  Filter  Bandwidth  (Ba): 

1 

1 

1 

1 

Hz 

Number  of  Pulses  Integrated  (n): 

10000 

10000 

10000 

10000 

Integration  Efficiency  (E-j(n)): 

0.03 

0.03 

0.03 

0.03 

Integration  Time  (T): 

1 

1 

1 

1 

sec 

Pulse  Compression  Ratio  (K): 

150:1 

150:1 

150:1 

150:1 

Target  Cross  Section  Area  (a): 

(1)  Transmissometers: 

125 

125 

125 

125 

m2 

(2)  Radars: 

10 

10 

10 

10 

m2 

Probability  of  Detection  (Pd): 

0.995 

0.995 

0.995 

0.995 

Probability  of  False  Alarm  (Pfa) 

:  10- 

■2,  10-4. 

10-6,  10-8, 

or  10 

-10 

Specific  Attenuation  (a): 

0,  5, 

10,  or  20 

dB/km 

*  See  Section  II I. C  for  exact  meaning  of  P^  for  each  system 


specific  attenuation  of  either  0,  5,  10,  or  20  dB/km 


For  the  six  types  of  instrumentation  discussed  in  Section  II,  we 
see  these  system  parameters  are  typical  of  a  system  using  paraboloid 
antennas,  Impatt  diode  power  sources,  and  a  corner  reflector  for  the 
calibrated  target. 

The  numerical  value  for  total  system  losses  and  noise  figure  are 
typical  of  an  experimental  system.  Since  all  programs  are  written  in  an 
interactive  form,  it  is  possible  to  easily  modify  input  parameters  to 
represent  virtually  any  kind  of  transmission  system. 

A.  Graphical  Results  and  Discussion. 

In  this  section,  graphical  results  for  each  of  the  six  transmission 
systems  will  be  presented  and  examined.  At  this  point  attenuation  will 
be  considered  as  an  independent  parameter.  The  effect  of  rain  and  fog 
will  be  Included  in  Section  V.B  and  V.C. 

1.  One-way  Transmission  and  Square-Law  Detection. 

Figures  12  through  16,  displayed  below,  show  the  typical  perfor¬ 
mance  of  a  system  with  separated  transmitter  and  receiver.  Here,  the 
receiver  uses  a  square-law  detector  as  discussed  in  Section  II.B.l. 
Figures  12  through  15  show  processed  signal -to-noise  ratio  versus  range 
for  different  values  of  frequency  and  and  specific  attenuation.  Figure 
16  shows  the  variation  in  range  as  a  function  of  specific  attenuation, 
a.  Because  this  transmission  system  Is  over  a  one-way  propagation  path, 
no  target  Is  Involved.  Therefore,  rather  than  specifying  probability  of 
detection  and  probability  of  false  alarm  as  in  a  radar  system,  calcula¬ 
tions  of  range  as  a  function  of  specific  attenuation  are  presented  for 
a  predetermined  processed  signal -to-noise  ratio  of  ID  dB. 


From  the  curves  displaying  siqnal-to-noise  ratio  as  a  function  of 
range,  it  can  be  observed  that  range  decreases  significantly  as  the 
frequency  is  increased.  Also,  notice  as  specific  attenuation  increases, 
the  signal -to-noise  ratio  E/N  at  the  system  output  not  only  decreases, 
but  its  rate  of  decrease  with  respect  to  range  is  significantly 
increased. 

In  Figure  16,  we  see  the  effects  of  specific  attenuation  and  fre¬ 
quency  on  range  performance.  As  expected,  when  frequency  is  increased 
range  performance  will  deteriorate.  Also,  as  specific  attenuation  in¬ 
creases  the  range  is  seen  to  be  greatly  reduced.  Even  under  the  worst 
possible  weather  conditions,  a  range  of  between  three  and  four  kilome¬ 
ters  can  be  achieved  from  the  two-ended  system.  All  graphs  presented 
in  this  section  are  derived  from  computer  program  SLDET2E  shown  in 
Appendix  C. 

In  these  graphs,  the  systems  are  compared  on  the  basis  of  equal 
specific  attenuation  at  all  frequencies.  A  more  basic  comparison  is  on 
the  basis  of  identical  meteorological  conditions.  This  is  deferred  to 
the  next  two  sections,  8  and  C. 
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12.  Signal -to-noise  ratio  versus  Range  for  the  double- 
ended  system  with  square-law  detection,  o-  0  dB/km. 
For  unspecified  system  parameters,  see  Table  I,  p.  53 


SE  (DB) 


f=  35  GHz 
f=  94  GHz 
f=140  GHz 


f=240  GHz 


1000.  2000.  3000.  4000.  5000.  6000.  7  TOO.  800Q. 

RANGE  (METERS)  \ 


Figure  14.  Signal -to-nolse  ratio  versus  Range  for  the  double- 
ended  system  with  square-law  detection,  a-  10  dR/km. 
For  unspecified  system  parameters,  see  Table  I,  p.  53. 
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Figure  16.  Range  versus  Specific  Attenuation  for  the  double- 

ended  system  with  square-law  detection.  E/N*  10  dB. 
For  unspecified  system  parameters,  see  Table  I,  p.  53. 
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2.  Two-way  Transmission  and  Square-Law  Detection. 

Figures  17  through  23,  presented  below,  show  typical  performance 
of  a  single-ended  system  with  the  receiver  Identical  to  the  previous 
system.  Again,  the  receiver  uses  a  square-law  detector  as  discussed  in 
Section  II.  ft. 2.  Specific  attenuation  Is  treated  as  an  independent 
parameter.  The  specific  effect  of  rain  and  fog  is  treated  in  Sections 
B  and  C. 

Figures  17  through  20  show  the  signal -to-noise  ratio  versus  range 
for  various  values  of  frequency  and  specific  attenuation.  Upon  compar¬ 
ison  of  Figures  17  through  2D  with  Figures  12  through  IS  of  the  previous 
system,  we  can  observe  this  single-ended  system  is  more  sensitive  to 
effects  of  specific  attenuation.  This  is  true  because  doubling  the 
effective  length  by  placing  transmitting  and  receiving  antennas  near 
each  other  causes  transmitted  signal  to  be  attenuated  twice  as  much  as 
for  the  double-ended  system. 

,  Figures  21  and  22  show  the  probability  of  detecting  a  target  both 
as  a  function  of  processed  signal -to-noise  ratio  and  range.  For  both 
graphs,  we  are  assuming  predetermined  false  alarm  probabilities  of  10-2, 
10-4,  io-6,  io-8,  anfj  io-1Dj  Figure  22  for  a  transmitted  frequency  of 
94  GHz.  From  Figure  21,  note  the  processed  signal -to-noise  ratio  must 
have  a  value  of  a£  least  18  dR  for  almost  certain  detection  of  the  125 
m2  target.  Under  the  assumptions  of  Table  I,  this  requires  the  target 
to  be  no  more  than  about  1300  meters  from  the  transmitting  and  receiving 
antennas  for  94  GHz. 

Finally,  Figure  23  shows  the  effects  of  specific  attenuation  and 
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frequency  on  radar  range  performance.  For  a  0.995  probability  of  detec 
ting  the  target  and  10"®  probability  of  false  alarm,  the  computed 
processed  signal -to-noise  ratio  at  the  output  of  the  detector  Is  15.83 
dB. 


As  in  the  case  of  the  double-ended  system  discussed  in  the  last 
section,  performance  will  decrease  as  frequency  and  specific  attenuation 
Is  increased.  All  graphs  in  this  section  are  derived  from  the  computer 
program  SLDET  shown  in  Volume  II. 
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Figure  17.  Slgnal-to-nolse  ratio  versus  Range  for  the  single- 
ended  system  with  square-wave  modulation  and 
square-law  detection,  a-  0  dB/km,  o=  125  m2. 

For  unspecified  system  parameters,  see  Table  I,  p.  53. 
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Figure  18.  Signal -to-nolse  ratio  versus  Range  for  the  single- 
ended  system  with  square-wave  modulation  and  square- 
law  detection,  a*  5  dB/km  ,  o*  125  m2. 

For  unspecified  system  parameters,  see  Table  I,  p.  53 
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Figure  19.  Signal -to-nolse  ratio  versus  Range  for  the  single- 
ended  system  with  square-wave  modulation  and  square- 
law  detection.  a=  10  dB/km ,  <j=  125  m2. 

For  unspecified  system  parameters,  see  Table  I,  p.  53 
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Figure  20.  Signal -to-nolse  ratio  versus  Range  for  the  single- 
ended  system  with  square-wave  modulation  and  square- 
law  detection,  a*  20  dB/km,  o-  125  m2. 

For  unspecified  system  parameters,  see  Table  I,  p.  53. 
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Figure  22.  Probability  of  Detectfon  versus  Range  for  the  single- 
ended  system  with  square-wave  modulation  and  square- 
law  detection,  f*  94  GHz,  o«  20  dB/km,  a*  125  nr. 

For  unspecified  system  parameters,  see  Table  I,  p.  53 
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[  23.  Range  versus  Specific  Attenuation  for  the  single- 

ended  system  with  square-wave  modulation  and  square- 
law  detection.  0.995,  Pfa»10"°  (Computed 
processed  E/N=  15.83  d8.) 

For  unspecified  system  parameters,  see  Table  I,  p.  53. 
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3.  CW  Signal  Transmission  and  Quadrature  Recei ver. 

Figures  24  through  30  show  typical  performance  of  a  radar  system 
which  transmits  a  CW  signal  toward  a  target  with  the  transmitting  and 
receiving  antennas  in  close  proximity.  However,  as  shown  in  Section 
II.B.3,  we  use  a  quadrature  receiver  rather  than  a  square-law  device 
as  in  the  previous  systems. 

Figures  24  through  27  show  the  signal -to-noise  ratio  versus  range 
for  various  values  of  frequency  and  specific  attenuation.  From  these 
graphs,  even  with  20  dB/km  attenuation  at  94  GHz,  a  range  between  1 
and  2  kilometers  is  expected.  This  is  a  significant  improvement  over 
the  two-way  transmissometer  system  using  square-law  detection. 

Figures  28  and  29  show  the  probability  of  detecting  a  125  m^ 
target  both  as  a  function  of  processed  signal -to-noise  ratio  and  range. 
For  both  graphs,  we  are  assuming  false  alarm  probabilities  of  10“2, 

10"4,  10-6,  10~8,  and  10"*0;  Figure  29  assumes  a  transmitted  frequency 
of  94  GHz.  As  with  the  previous  system,  to  be  certain  of  detecting  the 
target  a  si gnal -to-noise  ratio  of  about  18  dB  is  required.  A  signal -to- 
noise  ratio  of  about  18  dB  would  necessitate  the  target  to  be  no  more 
than  about  1700  meters  from  transmitting  and  receiving  antennas  for  a 
near  certain  probability  of  detection. 

Finally,  Figure  30  shows  effects  of  specific  attenuation  and 
frequency  on  radar  range  performance.  For  a  0.995  probability  of 
detection  and  10~6  probability  of  false  alarm,  the  computed  processed 
signal -to-noise  ratio  at  the  output  of  the  detector  is  15.83  dB.  All 
graphs  presented  in  this  section  are  derived  from  the  computer  algorithm 
C00ET1  presented  in  Volume  II. 
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Figure  24.  Signal -to-nolse  ratio  versus  Range  for  the  single- 
ended  CW  system  with  quadrature  receiver,  a*  0 
dB/km,  o=  125  m2. 

For  unspecified  system  parameters,  see  Table  I,  p.  53. 
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Figure  25.  Slgnal-to-nolse  ratio  versus  Range  for  the  single- 
ended  CW  system  with  quadrature  receiver.  a=  5 
dB/km,  a=  125  m^. 

For  unspecified  system  parameters,  see  Table  I,  p.  53. 
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Figure  26. 


Slgnal-to-nolse  ratio  versus  Range  for  the  single- 
ended  CW  system  with  quadrature  receiver.  a»  10 
dB/km,  o=  125  m^. 

For  unspecified  system  parameters,  see  Table  I,  p.  53 
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4.  Pulsed  Source  with  Incoherent  Pulse  Integration. 

As  discussed  in  Section  II. B. 4,  this  system  can  represent  either  a 
pulsed  source  range  gated  radar  or  a  transmissometer  detecting  the 
reflection  from  a  calibrated  target  such  as  a  corner  reflector.  Typi¬ 
cal  performance  of  this  radar  system  is  shown  In  Figures  31  through  37 
below. 

Figures  31  through  34  display  the  signal -to-nolse  ratio  versus 
range  for  various  values  of  frequency  and  specific  attenuation.  In  all 
cases,  10,000  pulses  are  integrated  with  an  integration  efficiency  of 

3%.  Since  these  graphs  are  for  a  pulsed  source  range  gated  radar,  a 

2 

target  cross  section  area  of  125  m  is  used.  Upon  comparison  of 
Figures  37  through  41  with  corresponding  figures  for  the  previous 
transmission  systems,  we  see  this  system  to  have  shorter  range  as 
specific  attenuation  is  increased.  However,  under  the  20  dB/km  propaga¬ 
tion  conditions,  minimum  ranges  of  between  800  and  1800  meters  can  be 
expected. 

Figures  35  and  36  show  the  probability  of  detecting  a  125  m2 
target  as  a  function  of  processed  signal -to-nolse  ratio  and  radar  range. 
For  both  graphs,  we  are  again  using  prescribed  false  alarm  probabilities 
of  10*2,  10“4,  10-6^  i0"8t  and  10-1^;  in  Figure  36  a  94  GHz  transmitted 
frequency  Is  assumed.  As  in  the  cases  discussed  before,  it  is  necessary 
to  have  a  processed  signal -to-noise  ratio  of  at  least  18  dB  for  a  near 
certain  probability  of  detecting  the  calibrated  target.  This  requires 
the  target  to  be  no  more  than  about  1000  meters  from  the  transmitting 
and  receiving  antennas. 
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Figure  31.  Signal-to-noise  ratio  versus  Range  for  the  single- 
ended  pulsed  source  range  gated  radar  employing 
post-detection  Integration.  n=10,000,  a*  0  dB/km, 
and  o=  125  m*. 

For  unspecified  system  parameters,  see  Table  I,  p.  53 
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Figure  32.  Signal -to-noise  ratio  versus  Range  for  the  single- 
ended  pulsed  source  range  gated  radar  employing 
post-detection  Integration.  (  n=  10,000  ,a«  5  dB/km, 
and  o=  125  m2. 

For  unspecifleo  system  parameters,  see  Table  I,  p.  53. 
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Figure  33.  Signal -to-nolse  ratio  versus  Range  for  the  single- 
ended  pulsed  source  range  gated  radar  employing 
post-detection  Integration.  n=10,000,  a°  10  dB/km, 
and  o=  125  wr. 

For  unspecified  system  parameters,  see  Table  I,  p.  53. 
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Figure  34.  Signal-to-noise  ratio  versus  Range  for  the  single- 
ended  pulsed  source  range  gated  radar  employing 
post-detection  integration.  n=10,000  ,a*  20  dB/km, 
and  o=  125  m?. 

For  unspecified  system  parameters,  see  Table  I,  p.  53. 


lV. 


82 


850.  900 


950.  1000.  1050.  1  100.  1  150.  1200. 

RANGE  (METERS) 


| — i — r 

8(?0. 


'  i  i 


Figure  36.  Probability  of  Detection  versus  Range  for  the 
single-ended  pulsed  source  range  gated  radar, 
f-  94  GHz  *aB20  dB/km,  and  cr=  125  m2. 

For  unspecified  system  parameters,  see  Table  I,  p.  53. 
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Finally,  Figure  37  shows  the  effects  of  specific  attenuation  and 
frequency  on  radar  range  performance.  For  a  required  probability  of 
detection  of  0.995  and  probability  of  false  alarm  of  10“®,  a  processed 
signal -to-noise  ratio  of  17.79  dB  is  computed. 

As  with  the  transmission  systems  analyzed  earlier  in  this  section, 
both  frequency  and  specific  attenuation  increases  can  greatly  reduce 
radar  range.  All  graphs  presented  in  this  section  are  derived  from 
the  computer  program  BOXCAR  given  in  Volume  II. 

5.  Pulsed  Source  with  £  Quadrature  Recei ver. 

As  discussed  in  Section  II. B. 5,  this  system  is  a  true  pulsed 
radar  which  uses  a  range  gated  quadrature  receiver.  Typical  performance 
of  this  radar  system  is  shown  In  Figures  38  through  44  presented  below. 

Figures  38  through  41  display  the  signal -to-noise  ratio  versus 
range  for  various  values  of  frequency  and  specific  attenuation.  From 
these  graphs,  note  that  even  with  20  dB/km  attenuation,  a  range  of 
between  1000  and  2000  meters  can  be  expected.  This  Is  slightly  better 
than  the  two-way  propagation  system  using  square-law  detection,  however, 
not  nearly  as  good  as  the  CW  source  with  quadrature  receiver. 

Figures  42  and  43  show  the  probability  of  detecting  a  10  m2 
target  as  a  function  of  processed  signal -to-noise  ratio  and  range.  As 
with  the  earlier  systems,  we  are  again  assuming  false  alarm  probabili¬ 
ties  of  10"2,  10'4,  10“6,  10"8,  and  10”10;  the  transmitted  frequency  In 
Figure  43  is  94  GHz.  For  this  pulse  radar,  a  processed  signal -to-noise 
ratio  of  at  least  18  dB  is  required  for  a  near  100%  probability  of 
detecting  the  calibrated  target.  This  will  require  the  target  to  be 
no  more  than  about  1950  meters  from  the  transmitting  and  receiving 
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S1gnal-to-no1se  ratio  versus  Range  for  the  single- 
ended  pulsed  radar  system  with  quadrature  receiver, 
a*  0  dB/km,c*  10  m2. 

For  unspecified  system  parameters,  see  Table  I,  p.  53 
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Figure  41.  Slgnal-to-nolse  ratio  versus  Range  for  the  single- 
ended  pulsed  radar  system  with  quadrature  receiver, 
a*  20  dB/km,  a=  10  m2. 

For  unspecified  system  parameters,  see  Table  I,  p.  53. 
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Figure  43.  Probability  of  Detection  versus  Range  for  a  single- 
ended  pulsed  radar  system  with  quadrature  receiver, 
f*  94  GHz  ><*=  20  dB/km,  and  o=  10  m2. 

For  unspecified  system  parameters,  see  Table  I,  p.  53. 
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antennas. 


Finally,  Figure  44  shows  the  range  versus  specific  attenuation 
for  the  pulsed  radar  and  quadrature  receiver.  For  a  0.995  probability 
of  detection  and  10"®  probability  of  false  alarm,  the  computed  processed 
signal  -to-noise  ratio  at  the  output  of  the  detector  is  17.79  dR.  All 
graphs  presented  here  are  derived  from  the  computer  algorithm  COOETB 
presented  in  Volume  II. 

6.  Pulsed  Source  with  Linear  FM  Rectangular-Envelope  Pulse  Compression 
As  discussed  in  Section  II. R. 6,  this  system  is  a  true  pulse  radar 
which  uses  a  pulse  compression  filter  and  range  gated  quadrature  recei¬ 
ver.  Typical  performance  of  this  radar  system  is  shown  in  Figures  45 
through  51  shown  below  where  a  pulse  compression  ratio  of  150:1  and 
10  m2  target  cross  section  area  are  assumed. 

Figures  45  through  48  show  the  signal -to-noise  ratio  versus  range 
for  various  values  of  frequency  and  specific  attenuation.  From  the 
graphs  of  signal -to-noise  ratio  versus  range,  even  with  20  dR/km 
attenuation,  a  range  between  1800  and  2600  meters  can  he  expected.  This 
is  about  the  same  range  as  predicted  with  the  single-ended  CW  system 
with  quadrature  receiver  for  a  comparable  radar  cross  section.  However, 
this  range  is  not  as  good  as  for  the  double-ended  system. 

Figures  49  and  50  show  the  probability  of  detecting  a  10  m2  target 
as  a  function  of  processed  signal -to-noise  ratios  and  radar  range.  As 
with  previous  systems,  a  transmitted  frequency  of  94  GHz  is  assumed  in 
Figure  50,  as  well  as,  false  alarm  probabilities  of  10"2,  10"4,  10"®, 
10"®,  and  10"10  in  both  figures. 
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Figure  45.  Signal -to-nolse  ratio  versus  Range  for  the  single- 
ended  pulsed  radar  svstem  with  pulse  compression, 
a*  0  dB/km  ,  o=  10  m*. 

For  unspecified  system  parameters,  see  Table  I,  p. 
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Figure  46.  Signal -to-nolse  ratio  versus  Range  for  the  single- 
ended  pulsed  radar  system  with  pulse  compression. 
a=  5  dB/km  »  <»■  10  m2. 

For  unspecified  system  parameters,  see  Table  I,  p.  53. 
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Figure  47.  Signal -to-nolse  ratio  versus  Range  for  the  single- 
ended  pulsed  radar  system  with  pulse  compression, 
a*  10  dR/km,  10  m2. 

For  unspecified  system  parameters,  see  Table  I,  p.  53. 
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Figure  48. 
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Signal -to-noise  ratio  versus  Range  for  the  single- 
ended  pulsed  radar  system  with  pulse  compression. 
cp  20  dB/km  ,  0=  10  m2. 

For  unspecified  system  parameters,  see  Table  I,  p.  53. 
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Figure  50.  Probability  of  Detection  versus  Range  for  a  single- 
ended  pulsed  radar  system  with  pulse  compression. 
f=  94  GHz,  «»  20  dB/km  and  a-  10  m^. 

For  unspecified  system  parameters,  see  Table  I,  p.  53. 
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SPECIFIC  RTTENUAT I  ON  (DB/KM) 


Figure  51.  Range  versus  Specific  Attenuation  for  the  single- 
ended  pulsed  radar  system  with  pulse  compression. 

Pd=  0.995,  Pfa=10-6  (Computed  processed  E/N=  17.79 
dB.) 

For  unspecified  system  parameters,  see  Table  I,  p.  53. 


Finally,  Figure  51  shows  the  range  versus  specific  attenua¬ 
tion  for  the  pulsed  radar  with  pulse  compression.  Using  the  same 
probability  of  detection  and  probability  of  false  alarm  as  In  previous 
systems  requires  a  processed  signal -to-nolse  ratio  of  17.79  dB  at  the 
output  of  the  detector.  All  graphs  presented  In  this  section  are 
derived  from  the  computer  algorithm  COMPRESS  presented  In  Volume  II. 


6.  Effects  of  Rain  Attenuation 

As  discussed  in  Section  IV,  It  Is  possible  to  develop  a  computer 
algorithm  to  compute  the  specific  attenuation  for  a  specified  rain  rate. 
For  this  study,  an  exponential  rain  drop-size  distribution  given  by 
Equation  (97)  and  a  routine  to  determine  complex  index  of  refraction 
for  liquid  water  Is  used  (see  Appendix  A.)  The  rain  water  refractive 
index  is  a  function  of  temperature  and  frequency.  Typical  values  of 
the  complex  Index  of  refraction  for  water  at  millimeter  wavelengths  are 
shown  in  TABLE  II  below. 

A  computer  subroutine  for  calculating  specific  attenuation  can  be 
applied  to  plot  maximum  range  for  any  of  the  six  systems  discussed  in 
Section  II  as  a  function  of  rain  rate.  In  this  calculation,  the  atten¬ 
uation  is  a  function  of  temperature;  however,  the  effects  of  temperature 
on  radar  range  is  negligible.  In  the  discussion  presented  here,  graphs 
showing  range  as  a  function  of  rain  rate  will  be  given  for  each  of  the 
systems  of  Section  II  with  a  temperature  of  20°C  assumed. 

The  drop-size  distribution  is  used  as  a  subprogram.  By  changing 

this  subroutine,  we  can  calculate  extinction  due  to  rain  with  other 

expressions  or  empirical  data  for  drop-size  distributions.  This  is 

useful  because  smaller  drops  become  very  important  in  the  millimeter 

wave  region  so  that  the  effect  of  drop-size  distribution  may  be  impor- 

« 

tant.  In  this  thesis,  the  Marshall-Palmer  exponential  distribution  [2] 
is  used  exclusively.  Also,  it  should  be  noted  that  the  subroutine 
works  well  for  attenuation  through  fog  when  given  an  appropriate  drop- 
size  distribution. 


TABLE  II 


Calculated  Complex  Index  of  Refraction  of  Pure  Water 


arious 


nnEsraa 


emperatures 


(GHz) 


Temperature  (°C) 


Complex  Index  of  Refraction 


35 

0 

4.032-j2.450 

35 

10 

4.610-j2.673 

35 

20 

5.200-j2.797 

35 

30 

5.771-j2.799 

94 

0 

2.810-jl.379 

94 

10 

3.073— j 1.653 

94 

20 

3. 359- j 1.929 

94 

30 

3.674-J2.191 

140 

0 

2. 575— j 1.015 

140 

10 

2. 748- j 1.254 

140 

20 

2.929-J1.514 

140 

30 

3. 131 -jl. 784 

240 

0 

2.405-j0.635 

240 

10 

2. 496-jO. 810 

240 

20 

2.571-jl.018 

240 

30 

2. 645- j 1.257 
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Figure  54.  Range  versus  Rain  Rate  for  the  single-ended  CW 
system  with  a  Quadrature  receiver.  P(j=  0.995, 

Pfa«  10-6,  a.  125  m2,  and  T=  20°C. 

For  unspecified  system  parameters,  see  Table  I,  p.  53. 
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RRIN  RATE  (MM/HR) 


Figure  55.  Range  versus  Rain  Rate  for  the  single-ended  pulsed 
system  with  a  quadrature  receiver.  Pd=  0.995, 

Pfa-  10-6,  o=  10  mz,  and  T=  20  °C. 

For  unspecified  system  parameters,  see  Table  I,  p.  53. 
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Figure  57.  Range  versus  Rain  Rate  for  the  single-ended  pulsed 
source  range  gated  radar  employing  post -detect Ion 
Integration.  P<j*  0.995,  Pfa«  10-6,o*  125  m2,  and 
T=  20  °C. 

For  unspecified  system  parameters,  see  Table  I,  p.  53. 


Figure  52  shows  range  as  a  function  of  rain  rate  and  frequency  for 
the  double-ended  system  with  square-law  detection  of  Section  V.B.l. 
Figure  53  shows  range  versus  rain  rate  for  the  single-ended  system 
with  square-law  detection  discussed  in  Section  V.B.2.  Figure  54  displays 
range  versus  rain  rate  for  the  single-ended  CW  system  with  quadrature 
receiver  system  of  Section  V.B.3.  Figure  55  shows  range  as  a  function 
of  rain  rate  for  the  single-ended  pulsed  system  with  quadrature  receiver 
discussed  in  Section  V.B.5.  Figure  56  gives  range  as  a  function  of  rain 
rate  for  the  single-ended  pulsed  radar  with  pulse  compression  discussed 
in  Section  V.B.6.  Finally,  Figure  57  shows  range  as  a  function  of  rain 
rate  for  the  range  gated  radar  of  Section  V.B.4.  As  expected,  the 
double-ended  systems  are  less  sensitive  to  the  effects  of  rain  because 
the  transmission  path  is  only  half  as  long  as  with  the  single-ended 
systems. 

C.  Effects  of  Fog  Attenuation. 

Not  only  can  we  construct  a  model  to  predict  radar  range  as  a  func¬ 
tion  of  rain  rate,  we  can  also  predict  radar  range  as  a  function  of 
liquid  water  content  in  fog.  As  discussed  in  Section  IV,  a  fog  drop- 
size  distribution  given  by  Equation  (108)  and  a  routine  to  calculate 
complex  Index  of  refraction  for  liquid  water  is  used. 

Here,  a  computer  subroutine  was  implemented  to  calculate  specif¬ 
ic  attenuation  of  fog  in  terms  of  liquid  water  content.  This  attenua¬ 
tion  calculation  can  be  applied  to  any  of  the  five  single-ended  systems 
discussed  In  before.  However,  here  we  choose  to  use  the  CW  radar 
system  of  Section  II. B. 3  in  order  to  compare  with  results  presented  in 
part  B  of  this  section. 
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Figure  58.  Range  versus  Liquid  Water  Content  In  fog  for  the 
single-ended  CW  system  with  a  quadrature  receiver 
Pd=  0.995,  Pfa=  10-6,  0=  125  m2,  and  T=  20  °C. 

For  unspecified  system  parameters,  see  Table  I,  p. 


53. 
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The  fog  drop-size  distribution  Is  used  as  a  subprogram.  By  chan¬ 
ging  this  subroutine,  we  could  calculate  attenuation  due  to  fog  with 
other  expressions  or  empirical  data  for  drop-size  distributions. 

Figure  58  shows  radar  range  for  the  CW  system  with  quadrature  re¬ 
ceiver  versus  liquid  water  content  of  a  fog.  This  graph  shows  frequency 
as  the  parameter  because  temperature  of  water  In  the  propagation  path 
does  not  significantly  affect  range.  The  temperature  and  frequency 
dependence  of  the  complex  Index  of  refraction  of  liquid  water  Is  shown 
In  Table  II,  as  computed  by  the  subroutine. 

In  general,  notice  that  the  radar  range  greatly  decreases  for  In¬ 
creasing  liquid  water  content  and  frequencies  In  fog  In  an  analogous 
manner  as  range  decreases  for  Increasing  rain  rates. 

D.  Conclusions. 

Six  transmission  systems  were  examined  for  frequencies  of  35,  94, 
140,  and  240  GHz.  The  first  system  examined  was  for  propagation  over 
a  one-way  path.  The  remaining  five  Instrumentation  systems 
measure  reflected  signals  from  a  calibrated  target.  Each  of  the  single- 
single-ended  systems  employs  a  different  method  of  signal  detection. 

It  Is  found  that  the  best  range  values  under  severe  weather  condi¬ 
tions  can  be  acquired  with  the  double-ended  propagation  system.  This 
Is  true  because  the  path  over  which  the  signal  travels  Is  only  half 
the  length  of  the  five  radar  systems  and  under  severe  conditions  path 
attenuation  can  be  very  large.  As  far  as  the  single-ended  systems  are 
concerned,  the  system  with  a  CW  signal  and  quadrature  detector  (Section 
II. B. 3)  appears  to  give  the  best  range  performance  under  severe  weather 
conditions  at  millimeter  frequencies. 
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For  an  experimental  transmlssometer,  the  single-ended  systems 
have  an  advantage  in  that  the  range  performance  can  be  varied  easily  by 
simply  moving  the  calibrated  target.  Thus,  it  may  be  best  to  use  a 
radar  system  with  a  quadrature  receiver.  This  type  of  instrumentation 
Is  preferable  because  of  the  linear  nature  of  the  detector  which  works 
extremely  well  for  low  signal -to-nolse  ratios  at  the  input  to  the  detec¬ 
tor. 

Finally,  an  examination  of  the  effects  of  rain  and  fog  attenuation 
was  presented  in  this  chapter.  The  purpose  was  to  show  quantitatively 
how  range  decreases  as  a  function  of  Increasing  rain  rate  and  Increasing 
liquid  water  contents  of  fog  as  a  function  of  frequency  and  system  con¬ 
figuration. 
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APPENDIX  A 


MODEL  FOR  COMPLEX  INDEX  OF  REFRACTION 
To  calculate  the  attenuation  parameter  a,  it  is  necessary  to 
develop  an  accurate  model  for  the  complex  Index  of  refraction  of  liquid 
water  as  a  function  of  frequency  and  temperature.  An  empirical  model  of 
the  complex  refractive  index  of  Ice  and  liquid  water  has  been  performed 
by  Peter  S.  Ray  [12]  with  corrections  given  by  Falcone  [13].  In  general 
this  model  is  applicable  from  -20°C  to  0°C  for  Ice  and  from  -20°C  to  50 
°C  for  liquid  water.  Also,  the  model  Is  found  to  be  applicable  for 
wavelengths  of  2ym  to  several  thousand  kilometers  in  Ice  and  from  2  ym 
to  several  hundred  meters  for  liquid  water.  In  the  presentation  below, 
a  general  discussion  of  Ray's  model  for  both  ice  and  water  will  be  gi¬ 
ven.  However,  for  our  work  only  liquid  water  will  be  required. 

For  this  model,  consider  data  given  In  terms  of  the  complex  rela¬ 
tive  permittivity,  c=e'-j€'1  where  e‘  is  the  relative  dielectric  con¬ 
stant  and  e'1  represents  the  effects  of  loss.  These  are  related  to  the 
complex  Index  of  refraction  by 


e'=nr2-ni2 

(Al) 

e"=2nrni  . 

(A2) 

Solving  these  two  simultaneous  equations  gives 

nr  =  |e'  +  |(e')2  +  (e")2]l/2jl/2 

(A3) 

Pi  =  . 

(A4) 

In  these  equations,  or  represents  the  real  part  and  nj  represents  the 
Imaginary  part  of  the  complex  Index  of  refraction  (n*nr-jni).  For 
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water,  the  complex  permittivity  Is  a  function  of  wavelength  according  to 


the  Debye  model 


e"  =  (  e<-ei )  ( Xs/X)  (A6) 

l+(Xs/X)* 


where  Is  the  high  frequency  dielectric  constant,  es  the  static 


dielectric  constant,  and  xs  the  relaxation  wavelength.  The  temperature 

dependence  of  es  Is  determined  experimentally  [12]  as 

es  =  78.54C  1.0  -  (4.579xl0-3)(T-25)  +  (1.19x10-5) (T-25)2 

- (2. 8x10-8 )(T-25)3]  (A7) 


where  T  Is  given  In  degrees  Celsius. 

The  equations  used  In  the  calculation  for  centimeter  absorption  are 
an  extension  of  the  Debye  theory,  modified  by  a  frequency  Independent 
conductivity  <»=12.5664xl08  mhos/cm.  These  equations  are  given  as 


e' 


(A8) 


and 


e"  ■ 


l-oo 

1 +?t*s^  1  n  ( aow/^y+fifl’x  pn'^o"  T  +  DT.WSxTiriO 


(A9) 


where  og  Is  a  spread  parameter  In  the  experimental  measurement.  The 
parameters  e^,  ao,  and  xs  are  solely  a  function  of  temperature  and  are 
given  by 


£i  *  5.27137  +  0.0216474T  -  0.00131198T2  (A10) 

an  -  -16.8129  +  0.0609265  (All) 

T+Z7T” 


and 


Xs  *  (3.3836xl0-4) 
where  all  wavelengths  are  In  centimeters. 


(A12) 


I 


K 


k 


l 


i 


i 
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The  complex  index  of  refraction  is  given  by  Equations  (A3)  and  (A4) 
when  Equations  (A10),  (All),  and  (A12)  are  substituted  into  (A8)  and 
(A9).  In  general,  the  complex  refractive  index  will  then  be  of  the  form 


n  =  nr-jni 


(A13) 


For  ice.  Equations  (A8)  and  (A9)  are  still  applicable,  however,  the 
defining  equations  for  e-j ,  a q,  Xs  and  o  are  different.  Once  again, 
results  are  found  experimental ly  as 


=  3.168 


oq  »  0.288  +  0.0052T  +  (2.3x10“4)t2 


-1.25xl04 

T79S59TT+27TT 


o  =  1.26e 


*s 


9. 


990288xl0‘5e 


1.32xl04 

1.986'§rr+273) 


(A14) 

(A15) 

(A16) 

(A17) 


and 

es  =  203.168  +  2.5T  +0.15T2  .  (A18) 

Substitution  of  these  last  five  expressions  into  Equations  (A8)  and 
(A9),  then  into  (A3)  and  (A4)  gives  an  expression  for  the  complex 
refractive  index  of  ice. 

The  FORTRAN  computer  program  REFRAC  shown  in  Volume  II  permits  a 
calculation  of  the  complex  index  of  refraction  as  a  function  of  a  given 
temperature  and  variable  wavelength.  Figures  A1  through  A4  show  typical 
results  of  the  complex  refractive  index  of  liquid  and  ice  water  for 
various  temperatures.  The  scales  for  the  refractive  index  were  chosen 
to  make  it  easy  to  compare  plots  for  water  and  ice  directly  with  plots 
given  by  Ray's  article  [121. 
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APPENDIX  B 


PULSE  COMPRESSION  RADAR 


For  the  transmission  system  under  study,  the  transmitted  signal  Is 
of  the  form 

f(t)  =  Acos[$(t)]  -T<t<T  (Bl) 

7  7 


where  <|>(t)  Is  the  phase  angle  of  the  transmitted  frequency  defined  by 
Cook  [18]  as 


$(t)  =  t^t+yt*- 

7 


(B2) 


and  u  Is  a  constant  given  by  y=Au>/T=2irAf/T.  Therefore,  the  Instantane¬ 
ous  frequency  Is  a  linear  function 

w  *  djKt )  =  Wc+ut  (B3) 


We  wish  to  determine  a  relation  between  Input  power  and  output 
power  of  the  pulse  compression  filter  so  that  the  signal -to-nolse  ratio 
before  detection  can  be  defined.  Using  the  method  outlined  by  Cook, 
let's  consider  a  unity  magnitude  linear  FM  rectangular-envelope  pulse 
compression  (l.e.,  A=l)  and  assume  the  generalized  compression  filter 
transfer  function  is 

cu 

H(w,u)  =  e  .  (B4) 

Therefore,  the  pulse  spectrum  Is  given  by  the  Fourier  transform  of  f(t) 
That  Is, 

"  -jult 

F(u>)  =  /  f (t)e  dt  (B5) 


-fcA.unVj'.V.A  -•  T 
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dt 


so 


T/2  -jut 

F(u)  =  /  cos£uct+ytzJe 


or 


F(u)  =  1 

7 


T/2  j£(uc-u)t+yt2j  T/2  -j[(uc+u)t+yt2] 

-T/2  6  ^  +  -T/2 


dll. 


(B5) 


Then,  if  we  define 


and 


T/2  j  fuct+ut^l  -jut 
F+(u,uc,u)  -I  eL  7  J  e  dt 
-T/2 


T/2  -jpuct+yt2)  -jut 
fr-(a|»uc,y)  =  /  e  L  7  e  dt 


it  follows  that 


-T/2 


F(u)  =  l[F+(u,uc,y)  +  Mu.c^.y)] 


However,  notice  that 


F_(oJ>tuc.y)  =  F+fu.-Uc-y) 


SO 


F (u)  =  l[F+(u,uc,y)  +  F+(u,-Uc,-u)l 


(86) 


(87) 

(88) 

(89) 

(810) 


and 

G(u)  =  F(u)H(u,y)  =  l[H(u,y)F+(u,uc,y)  +  H(u,y)F+(u,-uc,-y)].  (Bll) 

7 

Therefore,  the  real  time- varying  output  function  becomes 
9l(t)  *  l^-"1[H(u,u)F+(u,uc,y)]  +  l^_1[H(u,y)F+(u,-uc,-y)] 

7  7 


V 


g(t)  =  igi(t)  +  ig2(t)  . 

7  7 


(B12) 


We  can  examine  both  gift)  and  gzft)  separately  in  order  to  derive  the 
output  function.  By  definition  of  the  Inverse  Fourier  Integral 

gift)  =  -1[H(u>,y)F+.(u,uc,y)] 

»  jut 

=  1  /  H(u,u)F+(u,uc,u)e 

7t t 


00  jfu^-u)2  T/2  jr(“r-‘*))ti+uti2l 
91(t)  -  1  J  J  e  1  7  J 

7V  e  -T/2 


jut 

dtj  e  du 


which  can  be  rearranged  as 


T/2  [7 

9l(t)=  1  /  e 

7V  -T/2 


j (~ Mt  1 2+ urt  1  -*-Ur2  -  (ur+utl-ut)2' 

J.7  2u  ?jj 


x  /  e 


ju£u-(uc+uti-ut)^ 


Now,  using  the  substitution 


U  =  U-ju^Mtl-Mt,; 
/2y 


(B13) 


(B14) 


T/2  j[f1+“Ctl+r  ‘  .  jt 


we  see 


gift)*  /7m  /  e 
^  -T/2 1_ 

But,  we  can  show 

«  ju2  *  jw 

/  e  du  =  /[cosu2+js1nu2]du  *  /i'e  T 


/  e  du)  dti.  (B15) 


so  that 


T/2  jfyti 
91 (t )s  /7u  /  e  L7 
TT  -T/2 


9l(t) 


■e-1 


jf  uct-ut2+wl  T/2  jptti 


7  TJ  /  e 

-T/2 


which  reduces  to 


9l(t)=  /2y[ 


s1nut(T/2)  e 

'"mm 


j[u)ct-pt2+TTl 

>  L  7  ttJ  . 


(816) 


Now,  let's  perform  a  similar  Integration  to  find  g2 (t ) .  Knowing 
9?(t)=  $  “^(u.yjF+fu.-Wc-u)] 


=  1  /  Hfu.ujF+fu.-uc.-uJe  du 

Tir  -» 


92(t) 


j(^-w)2  ^T/2  -jj^(toc+u))ti+uti2‘J  jut 


7?  -« 


which  can  be  rearranged  to 

T/2  -jrutT2+urti+o)r2-(<^+uti-ijt)2 
92(t)~  1  /  e  [?  7y  2u  . 


7?  -T/2 


^oo  -ju^u-(uc+uti-ut)2J 


du  dti 


Now,  using  the  substitution  given  in  Equation  (B14),  we  see 

T/2  -irtiti2+urti-H»v2-(uf.-t-uti-ut)2l  «  -ju2 
92(t )=  /Zv  J  e  L2  "Tu  2p  J  J  e  du  dti. 

2w  -T/2  _  -«•  J 

But  we  can  show 

00  -ju2  »  -jw 

/  e  du  =  /  [cosu2-js1nu2!l  du  =  /we  7 


(B17) 


(B18) 


(819) 
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so  that 


T/2 

92(t )=  ftv  f 

TH  -T/2 


-  j  I” i£t  1 2+UCt  1  -^^-((dp+utl-Mt)2! 

e  [2  Zv  2u  J/ife 


•J* 

1 


dtl 


or 


92(0  = 


-jfaJct-ut2+nl  T/2 

L  ?  7j  / 

-T/2 


jvtti 

e 


dti 


which  reduces  to 

- -jruct-vit2+ir] 

92  (t)  =  /2y.fT/2]2  slnut(T/2)  e  L  7  TJ  (820) 

1  *  pt(T/2) 

Upon  substitution  of  Equations  (816)  and  (820)  into  (812)  gives 
g(t)  =  j2uCT/2]^s1niitjT/2)  cos^c-pt2+irJ  (B21) 

for  the  time-varying  output  function  where  u=2wAf/T  is  the  swept-fre- 
quency  deviation.  Figure  8  of  this  thesis  shows  the  relationship 
between  the  input  and  output  waveforms  of  the  pulse  compression  filter. 

Since  we  previously  assumed  the  input  amplitude  was  unity,  the 
output-input  peak  power  ratio  is  found  by  squaring  the  amplitude  of  the 
output  pulse.  This  implies 

=  2u[T/2]2  =  2(2*Af/T)[T/2]2  *  TAf  .  (822) 

n  1  * 

The  time  function  of  the  compressed  pulse  given  by  Equation  (821) 

Is  recognized  as  having  the  familiar  slncx  envelope.  The  pulse  width 
is  t=l/Af  when  measured  4  dB  down  from  the  peak  amplitude.  The  spacing 
between  the  first  zeros  of  this  envelope  is  2/Af.  The  Instantaneous 
received  frequency  Is  a  swept  linearly  (as  was  the  transmitted  fre¬ 
quency),  and  the  relative  peak  amplitude  is  /TAf’  as  shown  In  Figure  8. 
Use  of  these  results  gives  an  expression  for  the  output-input  peak  power 
ralo  In  terms  of  the  pulse  compression  ratio  K, 
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Notice  that  the  shape  of  the  compressed  pulse  Is  of  the  form 
slncx  and  not  rectangular.  Because  of  the  high  sldelobe  levels,  this 
may  not  be  the  most  desirable  pulse  waveform  for  some  radar  applica¬ 
tions,  however,  it  should  be  sufficient  for  our  application  of  estima¬ 
ting  the  effects  of  rain  and  fog. 

The  radar  equation  Including  losses  before  any  signal  processing 
at  the  receiver  is  given  by  Equation  (13)  for  a  single-ended  transmis¬ 
sion  system.  Therefore,  using  Equation  (B23)  gives  an  expression  for 
signal  power  as  -2aR 

t rr 

P0  =  PtKGtGrLtLrqX2  10  .  (B24) 


The  pre-detection  signal -to-noise  ratio  for  the  radar  transmisso- 

meter  Is  -2aR 

TT 

S  =  PtKGtGrLtLroX2  10  .  (B25) 

If  we  consider  a  pulse  compression  radar  with  a  quadrature  receiver 

system  as  in  Section  II.  B.5,  then  the  output  signal -to-noise  ratio  can 

be  found  from  Equation  (58)  to  be  -2aR 

TT 

E  *  PfKG*GrULroX2TdB*flO  .  (B26) 

if  (i-mToBF  ■ n 

This  result  is  presented  in  Section  II.  B.6  for  the  pulsed  source  with 
linear  FM  rectangular-envelope  pulse  compression. 
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